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ABSTRACT 

I present an homogeneous analysis of the transit light curves of fourteen well-observed 
transiting extrasolar planets. The light curves are modelled using jktebop, random 
errors are measured using Monte Carlo simulations, and the effects of correlated noise 
are included using a residual-permutation algorithm. The importance of stellar limb 
darkening on the light curve solutions and parameter uncertainties is investigated 
using five different limb darkening laws and including different numbers of coefficients 
as fitted parameters. The linear limb darkening law cannot adequately fit the HST 
photometry of HD 209458, but the other four laws give very similar results to each 
other for all transit light curves. In most cases fixing the limb darkening coefficients 
at theoretically predicted values does not bias the results, but does cause the error 
estimates to be too small. The available theoretical limb darkening coefficients clearly 
disagree with empirical values measured from the HST light curves of HD 209458; 
limb darkening must be included as fitted parameters when analysing high-quality 
light curves. 

In most cases the results of my analysis agree with the values found by other authors, 
but the uncertainties I find can be significantly larger (by factors of up to three). 
Despite these greater uncertainty estimates, the analyses of sets of independent light 
curves for both HD 189733 and HD 209458 lead to results which do not agree with each 
other. This discrepancy is worst for the ratio of the radii (6.7a for HD 189733 and 3.7 
a for HD 209458), which depends primarily on the depth of the transit. It is therefore 
not due to the analysis method but is present in the light curves. These underlying 
systematic errors cannot be detected from the reduced data alone unless at least three 
independent light curves are available for an individual planetary system. 
The surface gravities of transiting extrasolar planets are known to be correlated with 
their orbital periods. New surface gravity values, calculated from the light curve results 
and the stellar spectroscopic orbits, show that this correlation is still present. New 
high-precision light curves are needed for HD 149026, OGLE-TR-10, OGLE-TR-56, 
OGLE-TR-132 and GJ436, and new radial velocity curves for the XO-1, WASP-1, 
WASP-2 and the OGLE planetary systems. 

Key words: stars: planetary systems — stars: binaries: eclipsing — stars: funda- 
mental parameters 



1 INTRODUCTION 

The discovery of extrasolar planets was one of the 
great scientific achiev ements of the twentieth century 
IjMavor fc Queloal 19951 ). made possi ble by technologica l ad- 
vances in stellar spectroscopy (e.g. iButler et al]|l996l ). At 
this point over 200 exoplanets have been discovered through 
radial velocity monitorin g of their parent stars , making this 
technique very successful l|Udrv fc Santos 200?]) • The limit a- 
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tion of the method is that it is not possible to directly mea- 
sure any planetary properties: for each system only lower 
limits on the mass and orbital separation of the planet can 
be obtained. 

The solution to this problem arrived five years later, 
with the first detection of the transit of an extraso- 
lar planet across the dis c of its parent sta r, HD 209458 
jCharbonneau et all l200d : iHenrv et ail 120001 ). Since then 
over thirty transiting extrasolar planetary systems (TEPs) 
have been discovered, the majority by wide-field photomet- 
ric variability surveys. Much more information is available 
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for TEPs because the light curve shape is sensitive to the 
sizes of the star and planet and the inclination of the orbit. 

From the study of the light curve of a TEP, the spec- 
troscopic orbit of the star, and one additional piece of infor- 
mation (often a constraint on the stellar mass using theoret- 
ical models), it is possible to precisely measure the masses, 
radii, surface gravities and densities of both the star and the 
planet. This powerful procedure then allows detailed investi- 
gations into the chemical composition of the planet and the 
physical processes which govern its properties and evolution. 

The characteristics of the known TEPs have been de- 
rived in this way by many authors, but using several analysis 
techniques and a wide variety of observational and theoret- 
ical constraints. Considering only the light curve analyses, 
various studies have used different limb darkening functions 
and in many cases the coefficients of these functions have 
been fixed at theoretically predicted values without any con- 
sideration of whether this approach might introduce system- 
atic errors into the results. It is important to pay attention 
to the effects of limb darkening when modelling the light 
curves of TEPs because the coefficients may be correlated 
with the other parameters and so affect the measurement of 
the physical properties of the system as a whole. In light of 
this, it is appropriate to concentrate on studying the known 
TEPs to obtain their properties as accurately and as homo- 
geneously as possible. 

In this work I present light curve analyses for all four- 
teen TEPs for which decent light curves are available. Each 
system is studied using the same methods, to ensure that 
the results are as homogeneous and reliable as possible. A 
requirement of this method is that the importance of effects 
such as limb darkening must be quantified so the best ap- 
proach can be identified and used for every system. 

To enable an homogeneous study of the light curves of 
TEPs I investigate the following questions: Is linear limb 
darkening an acceptable approximation? Is it possible to 
measure accurate limb darkening coefficients from the data? 
Are these coefficients consistent with theoretical predic- 
tions? How accurately can orbital inclinations be measured 
from individual light curves? How important are correlated 
errors in each light curve? 

The results of this analysis makes it possible to answer 
further questions: Do the results in this work agree with 
published values? Do independent light curves of the same 
TEP agree with each other? If not, what is the size of the 
systematic error which causes the disagreement? How well 
understood are the properties of the known TEPs? 

In addition to dealing with these questions, I identify 
which of the studied TEPs would benefit the most from fur- 
ther study, both photometric and spectroscopic. Finally, I 
measure the surface gravities of each planet directly from the 
light curve parameters and the velocity a mplitude of the par- 
ent sta r, usi ng the method intr o duced by ISouthworth et all 
l|2004d ') and ISouthworth et"aH (|2007bl ). In a second paper 
(in preparation) the results presented here are used to de- 
termine the physical properties of the planets and stars in 
these systems using additional constraints from theoretical 
models, paying close attention to random and systematic er- 
rors and to the interagreement of independent sets of stellar 
evolutionary models. 



2 MODELLING THE LIGHT CURVES OF 

TRANSITING EXTRASOLAR PLANETARY 
SYSTEMS 

A TEP is simply a special case of an eclipsing binary sys- 
tem. The study of eclipsing binary light curves has been 
part of stellar astrophysics for many years, and was orig i- 
nally put on a solid foundation by Russell (| 1912al : Il912bh . 
Russell's technique, sometim es referred to as "rectification" , 
was s uccessively refined fsee lRussell fc Merrill|[l952l : iKopall 
Il959h but still requires several approximations which are not 
appropriate for many systems. 

When computers became widely available the Russell- 
Merrill method was superseded by models based on nu- 
merical integration, including the Wi l son-D evinne y code 
ilWilson fc Devinnevlll97]J:IWilsonlll979l.ll993D wink ilWoodl 



Il97ll . Il973h ebop (lEtzeJl 19751: iPopper fc Etzell fl981l) and 
others (|Hilllll97 9: Hadrava 1990 j) , wh ose u se is wide s pread . 
More recently, iMandel fc Agoll |2002l) and ICimenezI (120061 ) 
have produced fully analytical models which are intended for 
the study of TEPs. This approach allows a very fast compu- 
tation, but requires approximations such as the assumption 
of spherical shapes. 

One feature of a TEP system is that the radii of the two 
components are very different, which can be taxing for mod- 
els based on numerical integration over the surfaces of the 
components. Those models which use Roche geometry (e.g. 
the Wilson-Devinney code) simulate the component surfaces 
using several thousands of points each. This means that 
they very accurately reproduce distorted surfaces, but can 
be crippled by numerical error when the planet is so small 
that it only covers a small fraction of the points on the stel- 
lar surface. However, the use of a fully analytic model, with 
its accompanying requirement of sphericity, can potentially 
lead to an unquantifiable systematic error in TEPs where 
the planet is distorted because it is very close t o the star. 

In this work I adopt the ebop model (jEtzell Il975l : 
IPopper fc EtzelHigSlMEtzellligSlI ). which simulates the com- 
ponents of an eclipsing binary system using biaxial ellip- 
soids. It can also be restricted to spherical objects, so sys- 
tematic effects arising from the assumption of a physical 
shape can be easily quantified, ebop uses numerical integra- 
tion in concentric annu li over the surface of each component 
l|Nelson fc Davislll972l ). an approach which requires only a 
modest amount of computation and has a much lower nu- 
merical error than other non- analytic codes (|EtzelHl975r ). 

I use the jktebop codJ], a h eavily modified version 
of ebop (|Southworth et all 2004al lbT) which incorporates 
the Levenberg-M arquardt optimisation algorithm (mrqmin; 
iPress et al.l 1 19921 . p. 678) , an improved treatment of limb 
darkening (LD), and extensive error analysis techniques. 
The ebop model has been show n to work well for TEPs b y 
a number of a uthor s, including ISouthworth et all (l2007bh 
IWilson et all (|2006l ). ICimenezI (20061 ) and IShporer et all 
I 2007h . 



1 JKTEBOP is available from 

http : //www. astro .keele . ac .uk/~jkt/codes .html 
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2.1 What does a transit light curve tell us? 

An accurate orbital period, Porb, and reference time of min- 
imum light, To, is available for all systems studied in this 
work. Apart from P or b and To, the light curve of a TEP de- 
pends mainly on three quantities, which can be expressed 
in several different ways. In this work I take the quantities 
to be the fractional radii of the two stars and the orbital 
inclination, i. The fractional radii are defined as^\ 



where Ra and Rh are the (absolute) stellar and planetary 
radii and a is the orbital semimajor axis. There are three 
main directly measurable quantities: the transit depth, its 
overall duration, and duration of totality. Given a light curve 
with a well-defined transit shape, r\, ?"b an d i can be mea- 
sured accurately from these three main quantities. 

It is common in the literature to state that ordinary 
light curves do not reach the accura cy required to lift the 
degeneracy between i and Ra (e.g. iRowe et al.l 120071 ). A 
transit light curve has three main observable properties and 
depends on three main physical properties, so such a state- 
ment is incorrect. The crucial advantage the study of TEPs 
has over that of eclipsing binaries is that the optical light 
produced by the planet can be assumed to be zero, thus re- 
ducing the number of light curve parameters from four to 
three. The results obtained in this and other works demon- 
strates that ground-based light curves can provide definitive 
results for TEPs. Note that (apart from P or b) a light curve 
holds no information about the overall scale of the system 
so is not directly dependent on the mass or absolute radius 
of either component. 



values adopted here are not otherwise used in the deter- 
mination of the properties of the TEPs. To obtain a high 
numerical accuracy I set the size of the integration annuli 
(parameter 'intring' in ebop) to be 1°. LD is also included 
using one of five different laws (see Section l2~51l . 

One potential problem in analysing TEP light curves is 
that any additional light contribution from a nearby star will 
decrease the transit depth and cause a systematic error in 
the measured fractional radii of the two components. There 
is unfortunately no way of detecting 'third light' from a tran- 
sit light curve alone: including it as a fitted parameter will 
render the solution thoroughly indeterminate. Thankfully, 
the possibility of a blend is one of the basic tests performed 
by transit search pro jects to ensure a cand idate TEP is not 
a false positive fe.g. IWeldrake et al.ll2007h . The spectra of 
all confirmed TEPs have also received intense scrutiny and 
none have been found to contain light contamination from a 
third body. Third light can therefore be ruled out at roughly 
the 5% level for most systems. This amount of third light is 
in most cases insufficient to cause an important systematic 
error in the measured fractional radii. 

One of the assumptions of the light curve models used 
here and elsewhere is that the component bodies can be 
represented by geometrical shapes with definite edges. This 
assumption does not necessarily correspond to the canonical 
theoretical approach of defining a stellar or planetary sur- 
face as the point at which the optica l depth is r = |. Thi s 
effect has been considered in detail bv lBurrows et al.li|2007l ). 
who find that for a highly irradiated low-gravity planet the 
difference in radius can be up to 5%. This effect must be 
borne in mind during any comparisons between observed 
and theoretically predicted planetary radii. 



2.2 Solution procedure 

Rather than fitting for ta and directly, jktebop fits for 
the sum of t he fractional radii, (va + rb), and the ratio of 
the radii, k l|Southworth et al]|2005l ). This is because these 
quantities are only very weakly correlated for a wide variety 
of light curve shapes, which improves solution convergence. 
In the analyses presented in this work I have always fitted 
for (rA + r b ), k, i and T . 

The period, P or b, is generally known with much more 
accuracy than any other parameter, so I have usually fixed 
it at a known value. In most cases the orbita l eccentric- 
ity, e, is set to zero (see lLaughlin et al. 1 120051 for a well- 
presented study of previous claims of a non-zero eccentric- 
ity for HD 209458) and the planet is assumed to produce no 
light (an excellent approximation at optical wavelengths). 
The mass ratio, which in ebop governs how aspherical the 
planet is, is fixed at an appropriate value. The mass ratio 
has a negligible effect on the solutions (changes by a factor 
of two are unimportant). I emphasise that the mass ratio 

2 Throughout this work I identify stellar parameters with a sub- 
scripted 'A' and planet parameters with a subscripted 'b' to con- 
form to IAU nomenclature. I have neglected the usual subscripts 
's' and 'p' as these can lead to confusion (V for TEPs refers to the 
star but for eclipsing binaries has been used to designate the sec- 
ondary component; 'p' likewise could mean 'planet' or conversely 
'primary' in the different communities). 



2.3 Treatment of the stellar limb darkening 

LD is an important effect in light curve analysis of TEPs 
as it has a small but significant effect on the entire light 
curv e shape. This is often included using the simple linear 
law l|Russel]||l912bh : 

/ M = J„[l - u(l - /i)] (2) 

where [i — cos 7 and 7 is the angle between the surface nor- 
mal and the line of sight. 7o is the intensity at the centre of 
the stellar disc, J M is the intensity at angle 7 and u is the LD 
coefficient (LDC). The linear LD law is known to be a poor 
approximation to t he predictions of model atmospheres (e.g. 
IVan Hamme3ll993l b There is also evidence that it performs 
less well than more complex laws, at least at optical blue and 
visual wavelengths, i n eclipsing binaries ([Southworth et al.l 
1200781 ) and in TEPs (|Southworth et al.ll2007bl ). 

A currently unanswered question in the study of TEPs 
is how important is LD, and does the use of different laws 
make much difference to the light curve results? I have there- 
fore implemented a number of other LD laws into jktebop 
and used them in the ligh t curve mod elling process. These 
include the quadratic law (|Kopallll950l ): 

/M = M1-Ml- M )- Uq (l- M ) 2 ] (3) 

where the linear LD C is u a and the nonlinear LDC is v q ; 
the square-root law (|Diaz-Cordoves fc Gimenezlll992l) : 

/ m = /o[1-u s (1-m)-« S (1-v / A')] (4) 
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and the logarithmic law l|Klinglesmith fc Sobieskllll970h : 

Ip, = ^o[l — ui (1 — fj.) — vifiln/i] (5) 



Jl993h. iDfaz-Cordoyes etafl (ll995T). Idaret 



IClaret fc Hauschildtl <|2003h and lciaret et al 



I have also used the cubic law (Ivan 't Veerill960l : see also 
iBarban et afl|2003h : 



Jo[l-u c (l- M )-u c (l-/i) 3 ] 



(6) 



but this is complicated by the unavailability of theoretically- 
predict ed LDCs. 

In ISouthworth et all (|2007al ) it was found that the 
square-root and quadratic LD laws gave significantly bet- 
ter results than the linear law for an eclipsing binary with 
partial eclipses. Attempts to derive both linear and nonlin- 
ear LDCs for each star were successful, but gave imprecise 
results bec ause the two LDCs for o ne object are strongly 
correlated (|Southworth et all l2007al . their Fig. 3) . A simi- 
lar situation was consistently found in the current work (see 
Fig-El- We can draw three conclusions from this: (1) nonlin- 
ear LD laws should be used; (2) in most cases it is sufficient 
to fix one LDC at a reasonable value and fit for only the sec- 
ond LDC; (3) using LD laws with th ree or more L DCs (e.g. 
the four-coefficient law introduced bv lClaretll2000f ) is unnec- 
essary for data of the level of accuracy currently achievable. 

When modelling the light curves in this work I have 
calculated solutions for all five LD laws given above. I have 
always tried to fit for two LDCs when using a nonlinear law. 
In those cases where this did not give reasonable results, I 
have fixed the nonlinear LDC to a reasonable (theoretically 
predicted) value and fitted for the linear LDC. In a few cases 
this has also delivered unphysical results, so both LDCs had 
to be fixed to theoretical values. When an LDC was fixed, 
its value was perturbed in the error analysis to ensure that 
this source of error was included in the results. In every 
case the applied perturbations followed a flat distribution 
of ±0.1 around the fixed value. Although some nonlinear 
LDCs were fixed to theoretically predicted values, in these 
cases the linear and nonlinear LDCs were so strongly cor- 
related that this introduce only an infinitesimal dependence 
on theoretical model atmospheres. For a few systems I have 
also presented, for illustration and comparison with other 
studies, solutions with both LDCs fixed and also not per- 
turbed. 

Other researchers have fitted linear combinations of 
LDCs, fo r example 2w q + v„ and w q — 2v q for the quadratic 
law (e.g. IWinn et al.ll2007q ). These combinations are less 
strongly correlated than u q and « q so their values can be 
found more precisely. Using these combinations can improve 
convergence when fitting observational data, but I find that 
they do not lead to more precise measurements for u q and 
n q . I have not used these combinations for two reasons: no 
problems with convergence have been noticed; and u and 
v are more straightforwardly comparable to model predic- 
tions. 

Theoretically predicted LDCs were obtained by bi- 
linear interpolation in the stellar T c s and surface grav- 
itj0, using the tabulations provided by IVan Hammd 



3 The code for obtaining theoretically predicted LDCs by 
bilinear interpolation is called JKTLD and is available from 
http : //www. astro .keele . ac . uk/~jkt/codes .html 



2.4 Error analysis 

The formal parameter uncertainties produced by mrqmin 
are calculated from the solution covariance matrix and are 
unrel i able in the presence of param e ter correlations llPopper 
19841 : iMaceroni fc Rucinskil Il997l ; ISouthworth fc Clausen 
2007h . so have been ignored in this work. To calculate the 
random errors on the light curve parameters I have used 
the Mo nte Carlo simulation algo rithm implemented in jk- 
tebop (|Southworth et al . 2004b c), which w as found to give 
reliable results by Southworth et all (|2005h . 

The Monte Carlo algorithm finds the best fit to the 
actual data and uses this to generate a model light curve 
with the same time sampling. Then Gaussian noise is added 
and the resulting synthetic light curve is fitted in the same 
way as the actual data. The latter step is performed many 
times, and the 1 a uncertainties in the measured light curve 
parameters are assessed from the distribution of the best fits 
to the simulated light curves. A vital part of this process is 
that substantial random perturbations are applied to the 
initial parameter values before fitting each simulated light 
curve. This allows a detailed exploration of the parameter 
space and of parameter correlations, and so gives similar 
results to the Markov Chain Monte Carlo techniques used by 
other researchers (e.g. lHolman et ai1l2006l ). For each analysis 
I have used 1000 Monte Carlo simulations, which is a good 
compromise between avoiding small-number statistics and 
requiring a large amount of CPU time. 

Whilst the random errors are straightforward to quan- 
tify, the presence of correlated 'red' noise can cause system- 
atic errors which the Monte Carlo algorithm will not prop- 
erly account for. This is a widespread problem in the study of 
TEPs because changes in airmass, telescope pointing and fo- 
cus, detector response and other effects can cause significant 
correlations between adjacent datapoints in a light curve. 
Problems can also be caused by the data reduction tech- 
nique. Unfortunately, there is no general way of completely 
removing these effects from a light c urve. An exce l lent ex - 
ample of systematic bias is given by iGillon et al.l (|2007cT ) , 
who demonstrate that the depth of an observed transit of 
OGLE- TR- 132 varies depending on exactly how the data re- 
duction is performed. This bias could not be detected from 
the reduced light curve alone. A detailed investig ation of how 
red no ise affects the analysis of TEPs is given bv lPont et al.l 
(|2006h 

To assess the effect of systematic errors on the un- 
certainties of the light curve parameter measurements, I 
have implemented into jkte bop a residual-per mutation 
(or "prayer bead") algorithm (I Jenkins et al.l [20021 ). In this 
method, the residuals around the best-fitting light curve 
solution are all shifted to the next datapoint. A new fit 
is found, then the residuals are shifted again, with those 
which drop off the end of the dataset being wrapped around 
to the start of the data. This process continues until the 
residuals have cycled back to where they originated, and re- 
sults in a distribution of fitted values for each parameter 
from which its uncertainty can be estimated. The residual- 
permutation algorithm gives reliable and easily interpretable 
results, but its statistical properties are not totally clear and 
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it is only sensitive to correlated noise on timescales different 
to the eclipse durati on. The method has been used before 
on several TEPs (e . g. feouchv et al.ll2005al : IPont et al.ll2005l : 
iGillon et al1l2007al ; Knutson et al.ll2007al ). 

For each light curve I have used the residual permuta- 
tion algorithm to assess the importance of systematic errors. 
In those cases where it gives larger parameter uncertainties 
than the Monte Carlo algorithm, which only considers Pois- 
son noise, I have adopted the residual permutation results, 
which are sensitive to both Poisson and correlated noise. 
The parameter error estimates presented in this work also 
include a contribution from the effects of using different LD 
laws, and in every case are lcr errors. When observational 
errors are available I have used these to set the size of the 
Gaussian noise level for the Monte Carlo analysis. In cases 
where this resulted in a reduced \ 2 (Xred) °f greater than 
unity I have multiplied the Monte Carlo errors by U x r ^ d 
to compensate for t he underestimated obser vational errors 
<|Bruntt et al.ll2006l ; ISouthworth et al.ll2007al ). 



3 RESULTS FOR EACH SYSTEM 

Below I present an analysis of the available light curves of 
each TEP for which good or definitive data are available. 
For each light curve a set of solutions are tabulated for each 
of the five LD laws and for fitting either one or two LDCs. 
As these tables are quite bulky they have been shifted to an 
Appendix. For each system there is also a table comparing 
the final light curve parameters with those from previous 
analyses. 

The TEPs are discussed in rough order of increasing 
complexity, starting with the TrES, XO and WASP and 
HAT systems. The TEPs discovered using photometry from 
the OGLEJ3 project follow; these are in general a bit more 
difficult to analyse because they are fainter and so have less 
good light curves. GJ436 is then discussed on the basis of 
infrared observations from the Spitzer space telescope; this 
analysis is more complicated as the orbit is eccentric and 
the light produced by the planet is not negligible. Finally 
I investigate HD 149026, HD 189733 and HD 209458. Whilst 
the very shallow transit of HD 149026 makes it difficult to 
study, the brightness of the last two objects means that the 
available observations of them are of remarkable quality. 



3.1 TrES-1 

TrES-1 was discov ered to be a TEP b y the Trans- Atlantic 
Exoplanet Survey. lAlonso et ail (|2004l ) presented discovery 
light curves, follow-up photometry, and high-precision radial 
velocities which showed it to have a radius similar to that of 
Jupiter but a significantly lower mass. An excellent z-band 
light curve of TrES-1 was obtained by IWinn et~aT1 (|2007d ) 
and cov ers three con s ecutive transits. The analysis of these 
data by IWinn et all (|2007cT ) used quadratic LD, with the 
LDCs encouraged towards theoretically predicted values by 
the inclusion of a mild penalty term in the figure of merit. 
TrES-1 has also been studied bv lSteffen fc Agoll (|2005l L who 



4 http : //www . astrouw . edu . pl/~ogle/ 




-0.04 -0.02 0.00 0.02 

Orbital phase 

Figu re 1. Phased light curve of the transits of TrES-1 from 
W inn et all l l2007d ) compared to the best fit found using JKTE- 
BOP and the quadratic LD law with both LDCs included as fitted 
parameters. The residuals of the fit are plotted at the bottom of 
the figure, offset from zero. 



found no evidence of a third planet from an a lysis o f pub- 
lished times of transit, and bv lKnutson et ail (|200 7bL who 
were not able to detect a secondary eclipse from a ground- 
based L-band light curve. 

I have modelled the IWinn et al.l i|2007ch li ght curve of 
TrES-1 (1149 datapoints), adopting P orb = 3.030065 d and 
fitting for the midpoint of the second of the three transits. 
A circul ar orbit was adopte d, along with a mass ratio of q = 
0.00085 (|Winn et al.ll2007cr ). Solutions have been calculated 
for three scenarios: (1) LDCs fixed at theoretical values; (2) 
including the linear LDC as a fitted parameter, and fixing 
the nonlinear LDC at a theoretical value [and perturbing it 
during the Monte Carlo analysis to ensure this approach is 
reflected in the uncertainties] ; and (3) fitting for both linear 
and nonlinear LDCs. The solutions are given in Table lATI 

TrES-1 is unusual among the systems studied in this 
work in that the nonlinearity of its LD is very small. The 
lowest residuals are, as expected, found for solutions where 
both LDCs are allowed to vary. Solutions with the linear 
LDC as a fitting parameter but the nonlinear LDC fixed 
are in good agreement with solutions with both LDCs fit- 
ted. However, the solutions where the LDCs are held fixed 
at theoretical values differ in two way: the parameters are 
slightly different (and in the case of rb the discrepancy is 
significant) and the error bars are lower. 

The solutions in which at least one LDC is a fitted 
parameter are in good agreement with each other. As the 
values are well-defined even when both LDCs are fitted, I 
have used these in calculating the final solutions. The solu- 
tions with both LDCs fixed at theoretical values have been 
rejected, because we unfortunately cannot trust theoreti- 
cal LDCs (see Section l3.14|l . For the linear law, the best- 
fitting iiA does not agree with theoretical predictions. The 
best light curve fit for the adopted solution type and the 
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Table 1. Final parameters of the fit to the IWinn et al 1 <2007ch light c urve of TrES - 1 from the jktebop analysis. These have been 
compared to the properties of TrES-1 found bv lWinn et al, and lAlonso et al. Quantities without quoted uncertainties 

were not given by those authors but have been calculated from other parameters which were. 



This work 



Winn et al. (2007c 1 Alonso et al. (2004 



Sum of the fractional radii rA±rb 

Ratio of the radii k 

Orbital inclination (deg.) i 

Fractional stellar radius r\ 

Fractional planetary radius r\> 



0.1097 ±0.0022 
0.1381 ±0.0014 

88.67 ±0.71 
0.0964 ±0.0018 



0.1088 
0.13686 ±0.00082 
> 88.4 
0.0957 ±0.0014 



0.01331 ±0.00035 0.01308 ± 0.00023 



0.130. 
88.5 



+0.009 
0.003 
-1.5 
-2.2 



0.01279 



quadratic LD law is shown in Fig. [T] along with the observa- 
tions and the residuals. 

The residual perturbation results are in good agreement 
with the Monte Carlo results, which indicates that system- 
atic errors are not important for these data. The final re- 
sults from this study are given in Table[T] an d compared 
to pub lished studies. The results presented by IWinn et al.l 
(2007c) are in good agreement, but are closer to the solu- 
tions I find when fixing both LDCs. This is expected, be- 
cause W inn et afl l)2007ch included a penalty term in their 
figure of merit to enc ourage the LDCs to be clo se to theoret- 
ical values. Both the IWinn et al l i|2007d ) and lAlonso et al.l 
(2004) quote smaller errorbars than found in the present 
study. This is not surprising given the inclusion of additional 
error sources in the current work (e.g. the effect of different 
LD laws and the inclusion of the LDCs as fitted parame- 
ters), and confirms that the Monte Carlo analysis algorithm 
in jktebop is reliable in the presence of parameter correla- 
tions. 



3.1.1 Parameter correlations 

Fig. [5] shows example plots of the output of the jktebop 
Monte Carlo analysis algorithm for TrES-1. The first four 
panels show the joint distributions of parameter values for 
Ta, f'b, k, (rA+^b) and i. There is a clear correlation be- 
tween r\ and rb and a strong correlation between ta and 
k. The latter correlation is expected because k is very well 
determined (as it primarily depends on the eclipse depth). 
The parameters (rA+f*b) and k are only weakly correlated, 
vindicating their use as parameters of the fit. 

The lower half of Fig.[5]presents the Monte Carlo output 
for the LDCs for the four different nonlinear LD laws. Strong 
correlations are seen for all four laws, which is expected be- 
cause the light curves do not contain enough information on 
LD to obtain well-determined values for two LDCs; this situ- 
ation is very simila r to that for the partially- eclipsing binary 
system /SAurigae l|Southworth et a l. 2007a). The strongest 
correlation is seen for the square-root LD law and the weak- 
est correlation for the cubic law, but these correlations have 
no clear effect on the other light curve parameters (Ta- 
ble lAl[) . This means that the cubic LDCs are the best- 
determined, so the cubic LD law should be used if the LDCs 
are considered to be interesting in their own right. 

Fig-E] shows the very weak dependence of the derived 
ta and rb on the LDCs for the square-root law. Very little 
correlation can be seen, which confirms that the values of 
the LDCs have very little effect on those for the fractional 
radii of the components. The square-root law was used for 
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Figure 2. Scatter plots showing some results of the Monte Carlo 
analysis used to measure the uncertainties in the light curve pa- 
rameters. The upper four panels show results obtained using lin- 
ear LD: it can be seen that r\ and (rA+^b) are strongly correlated 
with each other in this representative case, whilst rb and i are also 
correlated. The lower four panels show the correlations between 
the two LDCs for each of the nonlinear LD laws: the correlation 
is the strongest for square-root and weakest for cubic. 
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Figure 3. Scatter plots similar to Fig. [2] but showing the weak- 
ness of the correlations between the LDCs and the fractional radii 
of the star and planet. The square-root LD law was used for the 
results in this plot as this law should show the strongest correla- 
tions. 



this comparison because it is expected to have the strongest 
correlations (see Fig.[5J). 



3.2 TrES-2 

TrES-2 was the second transiting extrasolar planetary sys- 
tem discovered by the Trans-Atlantic Exoplanet Survey 
lO'Donovan et al"1l2006l ) and is both larger and more mas- 
sive than Jupiter. Excellent z-band light curves of three 
transits were presente d and analysed by the TLC project 
dHolman et all l2007bT). who worked in close collaboration 
with Soz zetti et al. 1 12007m to obtain a highly accurate and 
consistent set of properties for the system. Its relatively low 
orbital inclination (i.e. high impact parameter) means that 
the fractional radii of the two components can be determined 
to a high accuracy. 

I h ave studied the light curve obtained bv lHolman et al] 
(|2007bh . which contains 1033 datapoints and is superior to 
those presented in the discovery paper (IO 'Donova n et al.l 
2006). I assumed a mass ratio of 0.0012 l)Holman et al.l 
2007bl : ISozzetti et al.ll2007f l and fitted for the orbital period 
as the Holman et al. data cover 13 orbital cycles. 

The results of the jktebop analysis are given in Ta- 
bleau The solutions where both LDCs are fitted return 
scattered results and very uncertain values for the coeffi- 
cients, so these have been rejected. The solutions with one 
fitted LDC are well-defined and in good agreement with each 
other (using the different LD laws), so have been adopted 
when calculating the final results. Solutions with both LDCs 
fixed give parameter values in good agreement with the 
adopted results, but the errorbars are substantially smaller 
(as also found for TrES-1). It seems that solutions with 
LDCs fixed at theoretical values not only can give biased 
results, but give error estimates which are artificially small. 



Table 2. Final parameters of the fit to the [Holman et al. (2007b) 
light curve of TrES-2 from the JKTEBOP analysis . These have been 
compa red to the properti e s of Tr ES-2 found bv lO'Donovan et al.l 
and lHolman et a l. (2007b). Quantities without quoted un- 
certainties were not given by those authors but have been calcu- 
lated from other parameters which were. 





This work 


O'Donovan 
et al. (20061 


Holman et 
al. (2007b 1 




0.1460 ±0.0042 


0.1429 


0.1475 


k 


0.1268 ±0.0032 


0.1273 


0.1253 ±0.0010 


in 


83.71 ±0.42 


83.90 ±0.22 


83.57±0.14 


rA 


0.1296 ±0.0038 


0.1267 


0.1311 ±0.0021 


rh 


0.01643 ±0.00046 
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0.0164 ±0.0004 
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Figu re 4. Phased li ght curve of the transits of TrES-2 from 
iHolman et ahl |2007bh compared to the best fit found using JK- 
TEBOP and the quadratic LD law with the linear LDC included 
as a fitted parameter. The residuals of the fit are plotted at the 
bottom of the figure, offset from zero. 



These solutions should only be used when there is no alter- 
native. 

For the uncertainties in the light curve parameters I 
have adopted those found using residual permutation, as 
they are slightly larger than those from the Monte Carlo re- 
sults. My final solution is given in Tablc[2]a nd compared t o 
those o f lO'Donovan et al l (|2006l ) and lHolman et all (|2007rl . 
As with TrES-1 I find parameter values in good agreement 
with those from other studies, but my uncertainties are 
somewhat larger (by factors of up to 3). When studying 
HD 189733 and HD 209458 (see below) I find that my so- 
lutions of different light curves do not agree to within the 
uncertainties I find. This indicates that my uncertainties are 
not overestimated, and so that the uncertainties given in pre- 
vious studies can be quite optimistic. Th is phenomenon is 
not pe culi ar to TrES-2 or the an alyses of lO'Donovan et al.l 
(|2006h and lHolman et all (|2007rJ ) but occurs for most of the 
TEPs studied in this work. 
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Figur e 5. P h ased light curve of the transits of XO-1 from 
iHolman et al.l d2006Tl compared to the best fit found using JK- 
TEBOP and the quadratic LD law with the linear LDC included 
as a fitted parameter. The residuals of the fit are plotted at the 
bottom of the figure, offset from zero. 



3.3 XO-1 

This was discovered to be a TEP by the XO project 
{McCullough et al"1 l2006) . A light curve was also obta i ned b y 
Super WASP- North and is presented bv lWilson et al.l {2006). 
The light curves from these two small-aperture large-field 
instruments are not d efinitive. Good light curves were ob- 
tained and studied by IHolman et all (2006) using a variety 
of telescopes and passbands. 

In this work I have analysed the two best datasets, 
which come from the FLWO 1.2 m (z-band, two transits, 
821 datapoints) and the Palomar 1.5 m (i?-band, 209 data- 
points) telescopes. I have adopted P Qr b = 3.9 41534 d and a 
mass ratio of 0.001 (|McCullough et al.l [20061 '). The results 
are given m Table|M] for the FLWO data: the Palomar 
dataset gives substantially larger parameter uncertainties so 
has only been used as a reality check. The solutions with 
two fitted LDCs give the lowest scatter, and well-defined 
coefficients (except for the square-root law), so have been 
used to calculate the final results. As with TrES-1, the best- 
fitting ua for the linear law does not agree with theoretical 
predictions. The residual-permutation uncertainties are in 
good agreement with the Monte Carlo results, indicating 
that correlated noise is not important for these data. 

The quadratic-law best fit is shown in Fig. [5] and the 
final light curve parameters are given in Table[3] The agree- 
ment between the parameters found here and those of 
IHolman et alj {2006,) is excellent, but agai n my uncertainties 
are si gnific antly greater. The results of iMcCullough et all 
(2006) and lWilson et all (|2006T l also agree well, taking into 
account the lower quality of the light curves available to 
them. 
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Figu re 6. Phased l ight curves of th e transits of WASP-1 from 
IShporer et ail l|2007f ) (grey points) and lCharbonneau et al.l (|2007T ) 
(black points) compared to the best fits found using JKTEBOP (red 
line for the Shporer data and blue line for the Charbonneau data) 
and the quadratic LD law with the linear LDC included as a fitted 
parameter. The residuals of the fits are plotted at the bottom of 
the figures, offset from zero. 



3.4 WASP-1 

WASP-1 was presen ted as a candidate TEP by 
IChristian et alj d2006l). and its pla n etary nature was 
confirmed by ICollier Cameron et al.l (|2007l ). A decent 
/-band light curve, contain ing 583 datapo i nts co vering two 
transits, was presented by IShporer et alj (|2007|T^I. A good 
z-ban d light curve was obtained bv Icharbonneau et al.l 
(2007) and contains 657 datapoints. 

I have analysed both the Shporer and Charbonneau 
light curves here, adopting P or b = 2 .519961 d and a mass 
ratio of 0.00075 (jShporer et all 120071 ). The results for the 
Shporer data are given in Table|X4] fitting for two LDCs 
does not cause a significant improvement in the results so I 
have used the solutions for fitting one LDC. From this point 
on, solutions with fixed LDCs have not been calculated as 
they are less reliable than other approaches. This is because 
the imposition of fixed LDCs is an artificial constraint which 



5 The data obtained by IShporer et al. (2007) are available in 
PDF format from the electronic version of that paper. They were 
converted into machine-readable format for this work using the 
UNIX command /usr/bin/pdf totext -layout. 
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Table 3. Final parameters of the fit to the Ifiolman et al 1 d2007bl ) light curv e of XO-f from the JKTEBOP analysis. Th ese have been 
compa red to the properties of XO-1 found by iMcCullough et al, (2006) and I Wilson et al, | ll2006h , as well as those of iHolman et al,l 
(2006). Quantities without quoted uncertainties were not given by those authors but have been calculated from other parameters which 
were. 



This work McCullough et al. (20061 Wilson et al. (2006) Holman et al. ( 20061 

r A + r b 0.1003 ±0.0022 0.1080 0.1084 0.1000 

k 0.1317±0.0019 0.134 ± 0.004 0.138 ± 0.020 0.13102 ± 0.00064 

i (°) 89.06±0.84 87.7 ± 1.2 88.9 ± 1.0 m.3lt° f s 

r A 0.0886 ±0.0019 0.0953 0.0953 0.0884 

r b 0.01166 ±0.00035 0.01273 0.01312 0.01159 



Table 4. Final parameters of the fit to the IShporer et al.l (|2007| ) and Icharbonne au et al. I (120071) li ght curves of WASP-1 from the 
JKTEB OP analysis. These have been compared to the properties of WASP-1 found by IShporer et all (120071) and ICharbonneau et al.l 
(2007). Quantities without quoted uncertainties were not given by those authors but have been calculated from other parameters which 
were. 



This work 
(Shporer data) 



This work (final result) 
(Charbonneau data) 



Shporer et al. (20071 Charbonneau et al. (2007 



r\ + r b 
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in 

n, 



0.200Z 



H1.015 
-0.013 

1033+ 0016 
u.iuoo_ Q 001g 

84.0 to 90.0 
al82 +o.on 



0.0187 



0.013 
+0.0013 
-0.0016 



0.1918 



+0.0101 
0.0064 



1045+ 0015 

u,lu * J -0 .0015 

86.5 to 90.0 



0.1737 
0.0182 



+0.0057 
-0.0089 
+0.0007 



0.1909 
0.1015 ± 0.0013 

89.8 ± 1.9 
0.1733 ± 0.0091 

0.0176 



0.1019 ± 0.0009 
86.1 to 90.0 



can both bias the best fit and decrease the derived parame- 
ter uncertainties. 

The residual-permutation algorithm gives uncertainties 
which are substantially greater than the Monte Carlo uncer- 
tainties: this problem seems to arise because the WASP-1 
has a high orbital inclination (low impact parameter) which 
is therefore poorly constrained. This difficulty also occurs 
when analysing the Charbonneau light curve fTable lA3|) . 
and means that i is only constrained to lie between 86° and 
90°. 

For the final results I adopt my solution of the Char- 
bonneau light curve with parameter uncertainties from the 
residual-permutation algorithm. The high-inclination nature 
of this system means its photometric parameters (particu- 
larly i) are not very well defined despite the availability of 
high-precision light curves. More extensive high-quality pho- 
tometric observations are needed for WASP-1. My solution 
of the Charbonneau light curve is in reasonable agreement 
with the solut ion of the Shporer lig ht cu rve and with the re- 
sults given bv lShporer et al.l <|2007l i and lCharb onneau et al.l 
l|2007l ). although the published results have lower uncer- 
tainty estimates and slightly smaller ratios of the radii. My 
solutions are compared to the observations in Fig. [6] 

3.5 WASP-2 

WASP - 2 was presented as a TEP candidate bv lStreet et alj 
d2007l) and confirmati o n of its nature was presen t ed by 
ICollier Cameron et all (|2007n . ICharbonneau et al.1 (<2007n 
obtained a good z-band light curve (426 datapoints) which 
I have analysed here. I ado pted P or b = 2.152226 d from 
ICollier Cameron et al.l (|2007t ) and a mass ratio of 0.001. In 
contrast to WASP-1, this system has a lower inclination 
(i = 84.8°) and therefore much better-defined fractional 



Table 5. Final parameters of the fit to the Charbonnea u et al.l 
l|2007r t light curve of WASP-2 from the JKTEBOP analysis. These 
have been compared to the properties of WASP-2 found by 
ICharb onneau et al 1 d2007l). I hav e not included the parameters 
from lColli^r*^amerori et al. I J2007t) as they are very uncertain due 
to the larger observational errors in the SuperWASP discovery 
light curve. 



This work Charbonneau et al. (2007) 

r A +r h 0.1409 ±0.0067 

k 0.1313 ±0.0017 0.1309 ±0.0015 

i (°) 84.83±0.53 84.74±0.39 

r A 0.1245 ±0.0058 

r b 0.01635 ±0.00093 



radii. This comparison makes it clear that it is much eas- 
ier to measure the radii of TEPs which do not have central 
transits. 

Table lA"6l contains the results of the jktebop analysis 
and demonstrates that allowing two LDCs to be fitted in- 
stead of one does not improve the quality of the fit at all 
and results in very small nonlinear LDCs. The residual- 
permutation method gives very similar uncertainties to the 
Monte Carlo analysis. For the final results I therefore adopt 
the solutions with one LDC fitted and with Monte Carlo un- 
certainties (Table Q - these results are i n reas onable agree- 
ment with those of ICharbonneau et al.l (|2007T ). The best fit 
is plotted in Fig. [7] 

3.6 HAT-P-1 

This ob ject was found to be a TEP system bvl Bakos et al. 
I|2007al ) from data taken by the HAT project (jBakos et al. 
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Figur e 7. Phased light curve of the transit of WASP-2 from 
ICharbonneau et al.l fl2007l) compared to the best fit found using 
JKTEBOP and the quadratic LD law with the linear LDC included 
as a fitted parameter. The residuals of the fit are plotted at the 
bottom of the figure, offset from zero. 



I2002I , |2004 ). Excellent light curves from the FLWO 1.2 m 
(z-band, 1709 observations), Lick 1.0m (Z-band, 1475 dat- 
apoints) an d Wise 1.0 m telescop es have been presented and 
analysed bv lWinn et~ai] |2007ei ). 

I have s tudied the FLWO and Lick observations from 
W inn et all (l2007ej) in this work, adopting a mass ratio of 
0.00045 IjBakos et al.ll2007al ) and including P or b as one of the 
fitted parameters. For both datasets the solutions with two 
fitted LDCs gave no significant improvement over those with 
one fitted LDC (TablelATl and Table lA"8|) , so the latter were 
used to calculate the light curve results. The residual-shift 
algorithm uncertainties were adopted as these are larger 
than the Monte Carlo method. The FLWO and Lick results 
are in excellent agreement (Table[6] Fig. [8} so they were com- 
bined using a weighted mean to obtain the final results. The 
final result is in good agreement with published analyses, 
and the photometric properties of HAT-P-1 are very well 
determined. 



3.7 OGLE-TR-10 

The transitin g natu re of this object was discovered by 
Udals ki et al.l (12002a) i n the c ourse of the O G LE-HI project. 
Konacki et alj l|2003bl ) and iBouchv et al.l (l2005aT ) subse- 
quently found a small radial velo city motion, and its e xo- 
planetary status was confirmed bv lKonacki et al l (120051 ). It 
was originally thought to be abnormally large for its radius; 
improved photometry has removed this anomaly but it re- 
mains one of the most bloated of the known TEPs. 

Dedicated pho tometric studies of O GLE- TR-10 have 
been p resented by IHolman et all (|2007al ) and iPont et al.l 



been p r 

(|2007cT ). The former study is based on an /-band light curve 
obtained using the Magellan telescope, which fully covers 
two transits separated by almost one year. The latter study 
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Figu re 8. Phased light curves of the transit of HAT-P-1 from 
W inn et al . (2007c) compared to the best fits found using JK- 
TEBOP and the quadratic LD law. The FLWO z-band data are 
plotted at the top and the Lick Z-band light curve is offset by 
—0.02. The residuals of the fit are plotted at the bottom of the 
figure, offset from zero. 



includes V- and P-band data covering only part of a single 
transit, observed with the VLT. 

The Magellan and VLT observations do n ot agree on 
the dep th of the transit. This is attributed by IPont et aD 
(2007c) to system atic errors ari s ing fro m the use of im- 
age subtraction bv lHolman et al.l (1200781 ). The OGLE data 
indicate a transit depth of ~1.9% (/-band), the Magellan 
data have a transit depth of 1.1% (/-band), and the VLT 
light curve transit depth is 1.4% (V and P). A version 
of the OGLE light curve, detrended using the algorithm 
of iKruszewski fc Semeniukl d2003|). gives a transi t depth in 

agreement with the VLT data (IPont et al.ll2007aV 

I have anal ysed both the M agellan I Holman et al.l 

l2007al ) and VLT (|Pont et al.ll2007d ) light curves for OGLE- 
TR-10, and separate solutions are given for both datasets. 
A mass ratio of 0.00056 and an in itial period of 3.101278 d 
was adopted (|Holman et al.| [2007a). For the Magellan data, 
I have included the period as a fitted parameter. For the 
VLT data the orbital period is a crucial quantity: because 
only half of one transit was observed, P or b is very highly 
correlated with To and with r\ and n>. Adopting a fixed 
P or b and To might therefore bias both the resulting va and 
rb and their error estimates. To avoid these problems I have 
taken six times of minimum light from t he literature (the 
fiv e times quoted bvlHolman et al . 2007a and the To given 
bv lKonacki etld1 l2005F and have included P or b and To as 
fitted parameters, jktebop treats both the times of mini- 
mu m and the light curve da tapoints as observational data 
fsee ISouthworth et al.ll2007al ). allowing an optimal and sta- 
tistically correct solution to be found. Uncertainties in all 
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Table 6. Final parameters of the fit to t he IWinn et al 1 <2007eh F LWO a nd Lick light curves of HAT-P-1 from the JKTEBOP analysis, 
compared to those found by IWinn et al, | ll2007eh and lBakos**etaT, (2007a). Quantities without quoted uncertainties were not given by 
those authors but have been calculated from other parameters which were. 



This work (FLWO) This work (Lick) This work (final) Bak os et al. f 2007a) Winn et al. (2007c) 

0.1026 ±0.0037 0.1056 ±0.0056 0.1035 ± 0.0031 0.1089 0.1044 

0.1125 ±0.0011 0.1102 ±0.0047 0.1124 ± 0.0010 0.122 0.11094 ± 0.00082 

86.37 ±0.31 86.10 ±0.38 86.26 ±0.24 85.9 ± 0.8 86.22 ±0.24 

0.0922 ±0.0033 0.0951 ± 0.0053 0.0930 ± 0.0028 0.0971 0.0940 ± 0.0028 

0.01037 ±0.00046 0.01048 ± 0.00047 0.01043 ± 0.00033 0.0118 0.0104 ± 0.0004 
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Figur e 9. Phased light c urves of t he transits of OGL E-TR-10 
from lHolman et al.l j2007al ) (top) and lPont et al.l ll2007d) (V mid- 
dle, R lower) compared to the best fits found using JKTEBOP and 
the quadratic LD law with the linear LDC included as a fitted 
parameter. The residuals of the fit are plotted at the bottom of 
the figure, offset from zero. 



observed quantities are included in the Monte Carlo analy- 
sis, so the error estimates also take into account the uncer- 
tainties in individual times of minimum light. 

For the Magellan data the best results were found by 
fitting for only one LDC (Table lA"9)l and the largest parame- 
ter uncertainties were obtained by the residual-shift method 
- the presence of significant correlated errors is unsurprising 
given that OGLE-TR-10 is in a crowded stellar field. The 
best results for the VLT also come from fitting only one 
LDC, but in these cases the residual-shift method indicates 
that correlated errors are not important. 



The solutions of the two VLT light curves are not in 
great agreement (Tables [A10I and IA 1 If) . but the differences 
are less than 1.7 a for each parameter. For the final solution I 
have taken a weighted mean of the VLT V and R light curve 
solutions, with uncertainties slightly increased to take into 
account the differences. This solution does not agree well 
with th e Magellan data, as was also found by IPont et al.l 
(|2007ch . The best fits are plotted in Fig. [9] where it can be 
seen that the two VLT light curves give noticeably different 
solutions, but both have clearly deeper transits than the 
Magellan data. The residuals of the best fit confirm that 
in each case the model provides a good representation of 
the observations, so any discrepancies are due to the data 
themselves. 

The agreement between our final solution (based on the 
VLT V and R data) and publishe d analyses is not good. 
The results of iKonacki et al.1 |2005t ) were based on the low- 
precision OGLE photometry, so it is not s urprising that 
these d isagree with later analyses. Those of iHolman et al.l 
i|2007al ) are in acceptable agreement with our o wn solution 
of the Magellan data, whereas the results from IPont et al.l 
i|2007ch are not in accord with our analysis. Whilst the solu- 
tion of the VLT data found here should be robust, OGLE- 
TR-10 would clearly benefit from further photometric study. 



3.8 OGLE-TR-56 

This object was the second TEP discovered, and the first us- 
ing the extensive photometric databa se built up by OGLE . 
Its eclipsing nature was discovered bv lUdalski et al. | |2 002b) 
and it was confirmed to be a planet by IKonacki et all 
(|2003ah from Keck radial velocity observations. Further ra- 
dial velocities, and r e-analyses of t he OG LE lig ht curve, have 
been p resented by iTorres et al.l {2004) and iBouchv et al.l 
|2005al ). 

Aside from the survey-quality photometry obtained by 
O GLE, the only light curves of OGLE-TR-56 were obtained 
bv lPont et al.l ( 2007d ) . Observations covering one full transit 
were taken in the V and R passbands using the VLT and 
FORS1 instrument. The two light curves have a very good 
photometric precision, particularly considering the apparent 
magnitude of OGLE-TR-56, but only a small number of 
datapoints. I have studied these two light curves, graciously 
provided by F . Pont, using P or b = 1.211909 d and a mass 
ratio of 0.001 (|Pont et al.ll2007d ). 

For both light curves, solutions with both LDCs free 
gave poor results (see Table lAT2l and Table lAl3)l . The final 
results have therefore been taken from the solutions with the 
linear LDC free and the nonlinear one fixed at theoretically 
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Table 7. Final parameters of the fit to the Magellan l|Holman et al.ll2007ah and VLT l|Pont et al.ll2007ch light curves of OGLE-TR-10 
from the JKTEBOP analysis. These have been compared to the properties of OGLE-TR-10 found in previous studies. Quantities without 
quoted uncertainties were not given by those authors but have been calculated from other parameters which were. 





Magellan 


VLT V 


VLT R 


VLT (this 


Konacki et 


Holman et 


Pont et al. 




(this work) 


(this work) 


(this work) 


work, final) 


al. (20051 


al. (2007a) 


(2007c) 




0.1468 ±0.0091 


0.188 ±0.012 


0.163 ±0.011 


0.175 ±0.010 


0.1260 






k 


0.1022 ±0.0024 


0.1174 ±0.0020 


0.1172 ±0.0008 


0.1172 ±0.0007 


0.127 ± 0.018 


0.0990 ±0.0021 


0.112 ±0.002 


in 


86.0±1.1 


83.00 ±0.94 


84.64 ±0.90 


83.87 ±0.69 


89.2 ± 2.0 


88.1 ±1.2 




rA 


0.1332 ±0.0081 


0.168 ±0.010 


0.146 ±0.010 


0.157 ±0.009 


0.1118 




0.183 ±0.018 


n> 


0.0136 ±0.0011 


0.0198 ±0.0014 


0.0171 ±0.0012 


0.0182 ±0.0011 


0.0142 







Table 8. Final parameters of the fit to the I Pont et al. I <2007ch light curves of OGLE-TR - 56 from the JKTEBOP analy si s. These have 
been co mpa red to the p r opertie s of the system found bv lKonacki et al. N2003al 1.l Torres et al. h2004h , lBouchv et al. 1 11200531) . iBouchv et all 
l|2005af) and lPont et al.l {2007c]). Quantities without quoted uncertainties were not given by those authors but have been calculated from 
other parameters which were. 





This work 


This work 


This work 


Konacki ct al. 


Torres et al. 


Bouchy et al. 


Pont et al. 




(V) 


(R) 


(final) 


(2003a') 


(20041 


(2005a1 


(2007c1 




0.275 ±0.038 


0.259 ±0.048 


0.269 ±0.030 


0.255 


0.242 






k 


0.1014 ±0.0061 


0.0907 ±0.0116 


0.0991 ±0.0054 


0.121 


0.124 


0.114 ±0.004 


0.101 ±0.002 


i c) 


78.8 ±3.0 


81.8±4.1 


79.8 ±2.4 


86 ±2 


81.0 ±2.2 


81 to 90 






0.250 ±0.034 


0.237 ±0.041 


0.245 ±0.026 


0.227 


0.215 




0.26 ±0.01 




0.0253 ±0.0041 


0.0215 ±0.0062 


0.0241 ±0.0034 


0.0276 


0.0268 







expected values (and perturbed during the Monte Carlo 
analysis) . For both datasets the residual-shift algorithm gave 
similar results to the Monte Carlo, indicating that correlated 
errors are not important. For the _R-band data the solution 
with the linear LD law is unphysical (ua > 1.0) - the LDC 
must be fixed when solving this light curve. In addition, the 
logarithmic and cubic LD law solutions are discrepant so 
have not been included in the final solution. 

The best fits to the two VLT light curves are shown 
in Fig. 1101 the two fits are not visually in good agreement. 
This can be seen in the final fits to each dataset given in Ta- 
ble[8] but the two solutions are actually consistent when the 
sizes of the errorbars are taken into account. I have therefore 
used the weighted average of the two for the final solution 
for OGLE-TR-56 (Tabled . A comparison with published 
results for this TEP shows reasonable agreement when the 
quality of the OGLE data are taken into account. This sys- 
tem would benefit from more extensive photometry. 



3.9 OGLE-TR-111 

This was al so found to be transiting by OGLE 

i|Udalski et all l2002al ) a nd has a low den sity. Its TEP na- 
ture was confirmed by IPont et alj ll2004h an d light curve 
studies have been published by IWinn et alj (l2007bl) and 
iMinniti et all (2007). The dataset obtained by I Winn et all 
|2007bl l is the best and is easily available so was chosen for 
analysis. It was obtained in the (presumably Cousins) / band 
and contains 386 observatio ns. I adopted P or b = 4.0144479 d 
and a mass ratio of 0.0007 (|Winn et al.ll2007bh (remember 
that the value of the mass ratio has extremely little influence 
on the light curve solutions). 

As for OGLE-TR-10 and OGLE-TR-56 no substantial 
gain was found in fitting for both LDCs (Table lA14|l . so the 



nonlinear LDC was fixed for the final solution (and per- 
turbed as usual during the error analysis). The residual-shift 
method gave slightly larger uncertainties, so correlated er- 
rors are important for these data. The final results are given 
in Tableland plo tted in Fig. |l l| and show good agreement 
with the study bv lWinn et alj (|2007bl l. 



3.10 OGLE-TR-132 

OGLE- TR- 132 was originally discovered as a candidate 



transiting syste m by lUdalski et al. | ([20031 1 and confirmed 

' l2004ll. , 



to be a TEP by IBouchv et all ||2004|1. A light curve from 
the VLT was presented by iMoutou et all (|2004|). The data 
were re - analysed using image deconvolution by Gillon et al.l 
(2007c), who found that the previous light curve was af- 
flicted by systematic errors from the data reduction. An ex- 
cellent discussion and examples of this problem are given by 
Gillon et al. 

I h ave solved the VLT light curve from iGillon et al.l 
( 2007c]) (274 observations), which have good precision but 
do not quite cover the full transit, finishing immediately be- 
fore the point of fourth contact. Unfortunately this lowers 
the accuracy of the results, and make correlated noise more 
important . I adopted a period of 1.689868 d and a mass ratio 
of 0.0009 (|Gillon et alj|2007d 1. 

Solutions with two fitted LDCs gave poorly defined co- 
efficients (Table lAl5[) so I have adopted the results of the so- 
lutions with the nonlinear LDC fixed. The values of ta and 
7"b are in good agreement between the different solution sets. 
The orbital inclination is not measured very accurately, and 
it is important to obtain a new light curve which fully covers 
the transit. The largest parameter uncertainties are found 
by the residual-shift method, and this is probably because 
the light curve does not fully cover the transit. Table flOl con- 
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Table 9. Final parameters of the fit to the IWinn et al ] ll2007bl) light cu rve of OGLE-TR-111 from the JKTEBOP analysis. These have 
been com pared to the propert ies of OGLE-TR-111 found by IPont et al, | | |2004) ,I Winn et al, | ll2007br i (from the same data as analysed 
here) and llVlinniti et all (|2007tl . 



This work 



Pont et al. 



Winn et al. 



Minmti ct al. 



TA + 
k 

in 

rA 



''b 



0.0952 ± 0.0044 
0.1316 ±0.0023 

88.11 ±0.66 
0.0842 ±0.0038 
0.01107 ±0.00067 



0.0943 
0.1209 
86.5 to 90.0 
0.0841 
0.01017 



0.132 ±0.002 
88.1 ±0.3 



0.09288 

1245+ 00050 
u.izio_ 0030 

87.0 to 90.0 

0.0926±°;°°5 

0.1028 




-0.04 



-0.02 



0.00 
Orbital phase 



0.02 



0.04 



Figur e 10. Pha s ed ligh t curve of the transits of OGLE-TR-56 
from IPont et al. (2007c) compared to the best fit found using 
JKTEBOP and the quadratic LD law with the linear LDC included 
as a fitted parameter. The upper (blue) curve is for the V band 
and the lower (red) for the R band data. The residuals of the fit 
are plotted at the bottom of the figure, offset from zero. 



tains the final parameters and the best-fit model is plotted 
in Fig.[H 

3.11 GJ 436 

The M dwarf system G J 436 is a unique object among the 
known TEPs as both the planet and star are by far the 
smallest and least massive; the system also has a currently 
unexplained eccentric orbit. This was known to be a plane- 
tary system from spectroscopic radial vel ocity measurements 
(IButler et al. 1 120041; iManess et alj \200H ) and transits were 
first detected by Gillon et al. ( 2007bl) . A good light curve of 



1.00 



0.99 



E 0.98 



0.97 




0.96 
-0.03 



-0.02 



-0.01 0.00 0.01 
Orbital phase 



0.02 



Figur e 11. Phas e d light curve of the transits of OGLE-TR-111 
from I Winn et al. (2007b) compared to the best fit found using 
JKTEBOP and the quadratic LD law with the linear LDC included 
as a fitted parameter. The residuals of the fit are plotted at the 
bottom of the figure, offset from zero. 



Table 10. Final parameters of the fit to the lGillon et alj l l2007ch 
light curve of OGLE-TR-132 from the JKTEBOP analysis. These 
have been compared to th e properties of OGLE -TR-132 found by 
iMoutou et al.l j2004h and I Gillon et al.l l|2007ch . Quantities with- 
out quoted uncertainties were not given by those authors but have 
been calculated from other parameters which were. 



This work 


Moutou et al. 


Gillon et al. 




(2004) 


(2007c1 



k 



0.231 ±0.023 
0.0937 ± 0.0026 
83.3 ±2.4 

0.211 ±0.020 
0.0198 ±0.0024 



0.0812 ±0.0017 0.0916 ±0.0014 



the tr ansit has been obta i ned using the Spitzer Space Tele- 
scope (|Gillon et al.ll2007al ; iDeming et"ldl2007al). Svitzer has 
also b e en used to study the secondary eclipse ( Deming et al.l 
l2007al; iDemorv et al.ll2007l ). 

Torres! l|2007l ) has measured the parameters of the sys- 
tem from its light and radial velocity curves, theoretical 
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Figu re 12. Phas e d light curve of the transit of OGLE-TR-132 
from [Gillon et all l l2007d ) compared to the best fit found using 
JKTEBOP and the quadratic LD law with the linear LDC included 
as a fitted parameter. The residuals of the fit are plotted at the 
bottom of the figure, offset from zero. 



Figure 13. Phased light curve of the transit of GJ436 from 
Spitzer compared to the best fit found using JKTEBOP and the 
quadratic LD law with both LDCs included as fitted parameters. 
The residuals of the fit are plotted at the bottom of the figure, 
offset from zero. 



stellar models and its measured absolute A"-band magni- 
tude and J — K colour index. He found that the star 
has a radius inflated by approximately 10% over that pre- 
dicted by the models, an anomaly which has previously 
been noticed for low-mass stars in eclipsi ng binary systems 
(|Ribasl l2003l . l200fj ; iLopez-MoraTesl 120071), single low-mas s 
field stars studied using interferometry ( Berger et al.| [2006). 
and th e low-mass secondary components of interacting bi- 
naries (|Littlefair et al.ll200a ). 

I hav e chosen to analyse the Spitzer light curve pre- 
sented by iGillon et ail l|2007al) . which contains 355 data- 
points covering one transit at a wavelength of 8/im. It is 
not possible to constrain the orbital eccentricity and lon- 
gitude of periastron, e and ui, from one transit so I have 
adopted the values e = 0.14 ± 0.01 jDemorv et all 120071 ) 
and uj = 351° ± 1.2° (|Maness et al.ll2007l). I have also taken 
Porb = 2.64385 d from | Gillon et al.l (|2007aT ) and a mass ratio 
of 0.00016 from lTorresI (|2007T ). An additional complication 
is tha t the brightnes s of the planet is not negligible at 8 /im. 
From lT5eming et al.l (l2007al . Figs. 1 and 2) I find a surface 
brightness ratio of approximately 7%. The surface bright- 
ness ratio has been fixed at this value; changes of a factor 
of two have a smaller effect on the radii than the final error 
bars. Similarly, changes of e and ui within their quoted error 
bars have a negligible effect on the final results. 

Attempts to fit for LD did not lead to realistic results 
(Table [A"T6)) . Whilst the photometry achieves a remarkable 
precision (0.7mmag per point), the observational scatter is 
still a significant fraction of the very shallow transit depth. 
The Spitzer data are taken at 8 /im where the effect of LD is 
minor, so for the final solutions it was acceptable to fix the 
LDCs to predicted values. The effects of correlated noise 
are not important. The final results are given in Table fTTl 
and shown in Fig. 1131 The agreement with literature values 
is good, but the solution does not lead to particularly ac- 



curate radii for either star or planet. Improved photometry 
is needed, and may be provided in the near future by the 
Transit Light Curve project fsee lWinnll2007l ). 



3.12 HD 149026 

HD 149026 b is a very anomalous objects, with one of the 
highest densities of the known TEPs. Its small size sug- 
gests that much of its mass lies in a rock/ice core, and 
HD 149026 has been the subject of several theoretical inves- 
tigations concerning its formation and chemical composition 
il Johnson et al.ll2006l; iBroeg fc WuchterllbOOTljMatsuo et all 
l2007l ; IZhou fc Linll2007l . e.g.). 

The TEP n ature of HD 149026 was discovered by 
ISato et all ||2005 ] ). and light curve s have also been pre- 
sented by (jCharbonneau et ai1l2006l ) from the FLWO 1.2 m 
telescopqj. However, the very shallow transit (depth 0.3%) 
means that its light curve parameters remain ill-defined de- 
spit e the reasonably good quality of the available photome- 
try. [WoffelaL| (|2007l ) have studied the Rossiter- McLaughlin 



6 Some of the Monte C arlo simulation results presented by 
Charbonncau ct al. | J2006I. their Fig. 2) are concentrated in a 
sharp peak at i « 90° whilst most of the results occupy a much 
broader distribution at i ~ 86°. A probably similar effect was 
seen whilst writing the Monte Carlo algorithm for JKTEBOP and 
is caused by the value of i becoming almost exactly 90° during the 
model fitting process. When this happens a numerical derivative 
will vanish (or an alternative solution will not be taken) if it is 
calculated symmetrically, because the symmetrical nature of the 
situation means that a solution with e.g. i = 89.9 is equivalent to 
i = 90.1. Thus the fit becomes stuck at i = 90° without physical 
justification. This e f fect ca n be avoided by not allowing i > 90°. 
ICharbonneau et alj | |2006| ) correctly identify their sharp peak at 
i RJ 90° as a numerical problem and do not use the relevant sim- 
ulations, so their results will be unaffected. 
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Table 11. Final parameters of the fit to the Spitzer light curve of GJ 436 from the JKTEBOP analysis, compared to those found by other 
authors. Quantities without quoted uncertainties were not given by those authors but have been calculated from other parameters which 
were. 



This work 



Gillon et al. 



Gillon et al. 



Doming et al. 



Torres f2007 



k 



(°) 



rb 



0.0791 ±0.0029 
0.08284 ±0.00090 
86.43 ±0.18 

0.0731 ±0.0027 
0.00605 ±0.00023 



0.082 ±0.005 
86.5 ±0.2 



85.90 



+0.19 
0.18 



0.0836 

0.0758 ± 0.0035 
0.00634 



0.0834 ±0.0007 

0.0750 ±0.0033 
0.00626 




2.430 



75 80 85 

Orbital inclination (degrees) 



90 



Figure 14. Variation of the roo t -mean -square of the residu- 
als of the fits to t he [ Sato et al] l|2005tl (top panel) and the 
ICharbonneau et alj ||2006|) FLWO a (middle) and r (bottom) light 
curves as a function of the orbital inclination. In each case only 
the quadratic LD law was used and the upper (blue) line is for 
fits where both LDCs were fixed and the lower (red) line is for 
fits where the linear LDC was fitted for. 



effect in the system and iBozorgnia et al] (|2006l ) have 
searched for the prese nce of lithium and potas sium in the 
planetary atmosphere. Irlarrington et ail (|2007h observed a 
secondary eclipse of HD 149026 with Spitzer and found a 
high brightness temperature at 8 /im of 2300 ± 200 K for the 
planet. 

I n this work I h ave solved the (b+y)/2 light curve 
from Sato et al. J2005I) an d the g- and r- band data from 
ICharbonneau et alj 1 20061 ). My initial solutions converged 
to values for the orbital inclination of either i « 81° or 
i « 90°. To investigate this, additional solutions were calcu- 
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0.99 




0.98 



0.97 



0.96 



-0.04 



-0.02 
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Figur e 15 . Phas ed light cu rves of the transits of HD 1 49026 from 
ISato et al] l|2005h (top) and lCharbonneau et all l|2006l ) (g middle, 
r lower) compared to the best fit found us ing JKTEBOP and the 
quadratic LD law. For the lSato et al] [1200511 data the linear LDC 
was fitted for whereas for the ICharbonneau et alj (120061 ) LCs all 
LDCs were fixed at theoretically expected values. The residuals 
of the fit are plotted at the bottom of the figure, offset from zero. 
All three datasets have been binned by a factor of 5 for display 
purposes. The differences between the best-fitting solutions for 
the three light curves are clear through the partial eclipse phases, 
and show that better data are needed. 



lated with i fixed to values between 90° and 70° in steps of 
0.2°. The quality of the fit for these additional solutions is 
shown in Fig |14l and the results for the three light curves are 
not in good agreement with each other o r with the analogous 
diagram presented bv lSato et al ] (|2005l . their Fig. 5). 

The solutions for the g light curve were very poorly- 
behaved: standard approaches tended to go to an ill-defined 
solution with i w 81°, but the additional solutions did not 
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Table 12. Final parameters of the fit to the lSato et alJ d2005l ) and lCharbonneau et al.l (|2006fl light curves of HD 149026 from the JKTEBOP 
analysis. These have been compared to the properties of HD 149026 found by these authors. Quantities without quoted uncertainties 
were not given by those authors but have been calculated from other parameters which were. 





This work 


This work 


This work 


This work 


Sato et al. 


Charbonncau ct al. 


Wolf et al. 




(Sato data) 


(FLWO g data) 


(FLWO r data) 


(final) 


(2005) 


(2006) 


(2007) 


k 


144+0-064 

U - L **_0.014 

047+ ' 005 
U - U£ " -0.002 


n 1 r:s+0-Oio 

U-15S_0.005 

n n^+o- 0023 
u.uiMci_ 0014 


147+0-011 
0.048_ 003 


147+0-014 

u - x * -0.008 


0.169 
0.0514 




0.507 ± 0.0015 


i (°) 


83.0 to 90.0 


86.0 to 90.0 


87.0 to 90.0 


86.0 to 90.0 


85.3 ± 1.0 


85.8±1;» 
0.0083 


86.1 ± 1.4 


'b 


n iQQ+o.06i 
0065+ 0039 

U.UUDO_ 00Qg 


u.iou_ 005 
0.0U8U_ 0004 


0.0068±g;ggg 


140+0- 012 
0.0068l«;°08 


0.161 
0.0082 





reproduce the corresponding fit parameters even when i was 
fixed at values near 81°. This indicates that the % 2 surface 
in parameter space has substantial structure, and these data 
have not been considered further. 

The solutions for the r and (b+y)/2 light curves are 
better-behaved and the resulting parameters are in rea- 
sonable agreement (even though the best fits are not; see 
Fig. ll5[) . The (b+y)/2 light curve is the best dataset, and 
it was possible to obtain meaningful values for the linear 
LDCs when the nonlinear LDCs were fixed. These solutions 
were therefore adopted as the final results for the (b+y)/2 
data. Analyses of the r light curve tended to give unphysical 
LDCs, so these had to be fixed at theoretically-predicted val- 
ues for the final solutions of these data. Attempts to fit for 
both LDCs for either light curve resulted in an impressive 
variety of physically impossible solutions. 

The nature of the transit of HD 149026 means that 
the inclination remains very poorly constrained. This un- 
certainty over which part of the stellar disc the planet 
eclipses causes the fractional radii of both components to 
be ill-defined, and the upper and lower confidence intervals 
for individual measurements are quite unbalanced. I have 
adopted the Charbonneau r outcome as the final results 
for HD 149026 (Table [L2j) . but have increased the error es- 
timates to reflect the difference between the r and (b+y)/2 
solutions. This choice results in more precise parameter val- 
ues, compared to the (b+y)/2 solution, but is dependent on 
theoretically-derived LDCs. However, any systematic error 
arising from the use of theoretical LDCs should be much 
smaller than the random errors of the current solution. 

The final solution presented here is not in good agree- 
men t with the i, ta and rb values found in pr e vious stud- 
ies jSato et all 120051 ; ICharbonneau et all 120061 ; IWolf et all 
120071 ). It corresponds to a higher (but still poorly-defined) 
i, which means smaller r\ and rb values are needed to give a 
transit of the observed duration. The substantially smaller 
rb found here indicates that the planet in the HD 149026 
system has an even larger density than previously thought. 
The only reasonable conclusion that can be drawn is that 
better light curves are needed before we can attain a precise 
understanding of the HD 149026 system. 

After the above analysis was performed, a nother study 
of HD 149026 became available. I Winn et all l|2007al ) have 
obtained new photometry and arrive at a high-i solution 
which is in good agreement with my own results. The frac- 
tional radii of the star and planet are still not well defined, 
so additional high-precision photometry remains a priority 
for observers. 



3.13 HD 189733 



HD 189733 was found to be a TEP bv lBouchv et al.l |2005br ) 
and is now one of the best-studied systems due to its bright- 
ness and relatively deep transit. It has a faint stellar com - 
panion separated by 11.2" on the sky jBakos et al. 2006a ). 
HD 189733 A has a detectable magnetic field ( Moutou et al.l 
120071 ) so is an interesting target for studying the activ- 
ity of planet host stars. A r otational period of 11 .953 ± 
0.009 d was determined by iHenrv fc Whinl (|2008T ) from 
photometry of its starspot-induced brig htness modulation. 
iHebrard fc Lecavelier Pes Etanga (|2006h found eclipses in 
its Hipparcos light curve, which has helped to refine P or b- 

Ground-based phot ometry of HD 189733 h as been pre- 
sented an d studied by Bakos et alJ (|2006bl ). IWinn et alJ 
(2006) and lWinn et al.l I (|2007dl ). The primary and secondary 
minima have also b e en studied from space using Spitzer 
(iDeming et al.1 I2006J ; iKnutson et~aH l2007dl ; iTinetti et ail 
2007) . In the ne ar future ne w light curves fro m HST 
jPont et al.ll2007bl ) and MOST dCroll et al.|[2007bl ) will be- 
come available. Finally, iBaines et al.l ( 20071 ) measured the 
angular diameter and thus linear radius of HD 189733 A, 
which means that it is possible to determine the physical 
properties of both star and planet using purely empirical 
methods (Paper II). 

I h ave solved the light curves presented by I Bakos et al.l 
l|2006bl 1 (620 datapoints), and the FLWO 1.2 m z-band data 
(1662 observation s) and Wise 1 .0 m 7-band light curve (344 
datapoints) from IWinn e t al. (20 06|). A slight eccentricity 
has no t been ruled out (|Deming et al.ll2006l ; IKnutson et al.l 
l2007dl ) but such a small value has a negligible effect on 
the parameters d erived in thi s work . I have taken P or b = 
2.2185733 d from I Winn et all (|2006r i and a mass ratio of 
0.0014 from I Winn et al.1 (|2006h . 

For each of the light curves it was difficult to extract 
reliable values for both LDCs so I have adopted the so- 
lutions with one LDC fixed (Tables lA"l8l |A"T1)1 and lAlO)) . 
The Monte Carlo and residual-permutation uncertainties 
were similar so correlated noise is not important for these 
datasets. The best-fitting light curves are compared with the 
data in Fig. 1161 

The solutions of the three light curves are inconsistent 
(Table fC!)) , due mainly to a mismatch of 6.7ct in the ratio of 
the radii, k. This results in discrepancies of 3.1 ct in rA and 
2.1 ct in rb. As k depends essentially on the eclipse depth, 
this problem is related to the data rather than to the anal- 
ysis. It is not very surprising that this discrepancy exists, 
because obtaining light curves of TEPs is an observation- 
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Table 13. Final parameters of the fits to the iBakos et aD j2006bl ^ (FLWO r) and IWinn et all d2007dl l (FLWO z and WISE I) light 
curves of HD 189733 from the JKTEBOP analysis. These have been compared to the properties of HD 189733 found by other workers. 
Quantities without quoted uncertainties were not given by those authors but have been calculated from other parameters which were. 





This work 


This work 


This work 


This work 


Bouchy et 


Bakos et 


Winn et 




(FLWO r data) 


(FLWO z data) 


(WISE I data) 


(final result) 


al. (2005b) 


al. f2006b1 


al. C2007d) 




0.1184 ±0.0022 


0.1297 ±0.0023 


0.1326 ±0.0083 


0.1287 ±0.0036 


0.1311 




0.1301 


k 


0.1542 ±0.0019 


0.1577 ±0.0010 


0.1499 ±0.0029 


0.1568 ±0.0024 


0.172 ±0.003 




0.1574 


i (°) 


86.50 ±0.59 


85.73 ±0.20 


85.38 ±0.68 


85.78 ±0.25 


85.3 ±0.1 


85.79 ±0.24 


85.76 ±0.29 


r A 


0.1026 ±0.0054 


0.1121 ±0.0019 


0.1153 ±0.0071 


0.1113 ±0.0031 


0.1118 




0.1124 ±0.0034 




0.0158 ±0.0012 


0.0177 ±0.0004 


0.0173 ±0.0012 


0.0175 ±0.0005 


0.0192 




0.0177 ±0.0007 
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Figur e 16. Phased ligh t curves of the transits of HD 189733 
from IBakos et aLI <2006bD (top) and IWinn et all l l2007dh (mid- 
dle: FLWO z; bottom: WISE /) compared to the best fit found 
using JKTEBOP and the quadratic LD law with the linear LDCs 
included as fitted parameters. The residuals of the fit are plotted 
at the bottom of the figure, offset from zero. 



ally demanding activity. An example of systematic problems 
with the ecli pse depth in OGLE- TR-132 was presented and 
discussed by iGillon et~al] (|2007ch . For HD 189733 there are 
three high-precision light curves available, so the effects of 
systematic errors can be estimated in detail. HD 189733 is 
therefore one of the best-understood of the TEPs. A rather 
unhappy implication of this discrepancy is that we can only 
fully understand the intrinsic properties of a TEP if three in- 
dependent light curves are available. It is possible that the 
presence of starspots contributes to the different k values 
found for the light curves; this could be investigated using 



th e excellent HST ligh t curve which has since been presented 
bv lPont et all (|2007bl ). 

For the final solution I have calculated the weighted 
mean for each parameter from the three light curves stud- 
ied here, and have increased the uncertainties to account for 
the discrepancies discussed above. Because of this, the re- 
sulting solution is very robust. More pieces of the jigsaw can 
be added to this study in the near future by including the 
forthcoming HST and MOST datasets. The study of sys- 
tematic biases in light curve parameters is a vital part of 
understanding the extent of our knowledge of the properties 
of TEPs. 



3.14 HD 209458 

In clear contrast to its position as the final ob- 
ject considered in this work, HD 209458 was the first 
known TEP, with disco very light curves publi shed by 
ICharbonneau et all l|2000t ) and iHenrv et all (feoOCh . It has 
since become the most heavily studied, and the observa- 
tional data available for this system constitutes an em- 
barrassment of riches . Transit light cur ves have been ob- 
tained by HST (twice: Brown et alj|200ll andlKnutson et al.l 
l2007d). MOST dRowe et al.ll2006l : ICroll et alj2007al ). Svitzer 



Richardson et al. 1 boo lTa^ from terminal facilities 



( Wittenmver et aD 2005 ). Secondary eclipses have been ob - 
tained by Spitzer ( Deming et a l. 2005; Knutson et al.| [2007a) 
but gro und-based attempts have been onl y partially suc- 
cessful (|Snellenl [20051 : iDeming et alj|2007bh . Spectroscopic 
observations of the Rossiter-McLaughlin effect through 
transit have added to the und e rstan d ing of HD 209458 
jQueloz et al.1 1200(1 IWinn et all 120051 : IWittenmver" et al.l 
120051 ) and no less than three works have obtained precise 
P Qr b measurements using Hipparcos data |S5derhielmlll999l : 
iRobichon fc Arenoullioool : ICastellano et al.ll20od ) 

Given such a large amount of data, I have decided 
to analyse both set s of H ST light curves (|Brown et al.l 
l200ll: iKnutson et all l2007d ) and the MOST light curve 
i Rowe et al.ll2006l ). Studying these datasets allows careful 
cross-checks of the success of the different solutions and the 
interagreement of the results, and so indicate whether or 
not systematic errors (arising particularly from red noise) 
are important. This is particularly relevant in light of the 
results for HD 1897 33 presented above. The HST data of 
iBrown et al.l (l200ll) have been p reviously solved using JK- 
TEBOP ( Southworth et all l2007bl ) but are revisited here in 
order to ensure homogeneity of the results. I have adopted 
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Figu re 17. Phased light curves of the transi t of HP 209458 from 
iBrown et all ll200ll l (upper) and lRowe et al.l ||2006|) (lower) com- 
pared to the best fits found using JKTEBOP and the quadratic LD 
law. Each consecutive set of 10 points of the MOST light curve 
have been combined to make the plot clearer. The HST light curve 
is unbinned. The residuals of the fit are plotted at the bottom of 
the figure, offset from zero. 



P orb = 3.524 7 4859 d and a mass ratio of 0.00056 from 
iKnutson et "all i|2007d ). 

The HST data obtained by IBrown et all (|200ll ) (556 
datapoints) give good values for both LDCs when the non- 
linear LD laws are used (Table lA"2~lj) . so these solutions have 
been adopted. The high quality and relatively blue pass- 
band of these data allow the linear LD law to be ruled out 
at a high significance level (greater than 99.99%). There is 
less to choose between the nonlinear laws, except that the 
quadratic law gives marginally the best fit. The low Poisson 
noise means that correlated errors dominate here, so I have 
adopted the uncertainties from the resid ual-shift algorithm . 

When analysing the MOST data |Rowe et all 120061 ) 
I rejected the datapoints which were distant from transit 
(Fig. ll7p as these add nothing to the solution except in- 
creased computational requirements. I also did not consider 
data from the second of the four transits as only part of this 
eclipse was observed. This left a total of 8772 datapoints. 
Solutions with two fitted LDCs were not great so I have 
adopted those with one fitted LDC. As with the HST data, 
correlated noise is important so I have adopted the uncer- 
tainties from the residual-shift algorithm. 



3.14.1 HST observations of \Knutson et al\ \200lA) 

The be autiful HST observations obtained bv lKnutson et ail 
l|2007d ) allow the properties of the HD 209458 system to be 
determined to a very high precision. These data comprise 
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Figur e 18. Phased lig ht curves of the transit of HD 209458 from 
IKnutson et alj d2007ch compared to the best fits found using JK- 
TEBOP and the quadratic LD law (both LDCs fitted). The data 
are plotted with the bluest passband at the top and the reddest 
at the bottom, with offsets to make the plot clear. The residuals 
of the fits follow the same sequence and are plotted at the bottom 
of the figure with offsets from zero. 
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two separate sets of observations obtained with a blue (504 
epochs) and a red grating (548 epochs) using HST/STIS. 
Each set has been split into five different passbands, result- 
ing in ten separate light curves evenly distributed between 
320 nm and 980 nm. The variation of LD strength with wave- 
length is clear and can be a strong test of the predictions of 
stellar model atmospheres. 

The ten light curves have been solved separated using 
linear and quadratic LD (Tables IA23I and IA24|I . and the 
best fits and residuals are plotted in Fig. [18] The light curves 
have a low photon noise so the residual-shift algorithm gives 
the largest u ncerta inties. As with the HST observations of 
I Brown et al.l (|200l[ ). the linear LD law gives a significantly 
worse fit than the quadratic LD law for the bluer wavelength 
bands. In particular, the referee has correctly pointed out 
that, for the linear-law case, the orbital inclination values are 
correlated with the wavelength of observation. This is clearly 
unphysical and arises from the inadequacy of the linear law 
being more important at shorter wavelengths where LD is 
stronger. No correlation is seen for the quadratic-law case. 

The variation with wavelength of the linear LDC is com- 
pare d in Fig . ll9lto predictions for the SDSS ugriz passbands 
from IClaretl (|2004f ) . The disagreement between observation 
and theory is clear. However, detailed conclusions should not 
be drawn from Fig. [T5] for two reasons: the linear law is not 
a good representation of the LD of stellar discs, and this 
is more important at bluer wavelengths; and the plotted 
LDCs are not monochromatic but refer to different pass- 
bands. I find that the dependence of the light curve solu- 
tions on wavelength, noted in the previous paragraph, does 
not have a significant effect on the conclusions drawn from 
Fig.rjl 

In Fig. [20] the observed and predicted LDCs are com- 
pared for the quadratic law. The observed LDCs again do 
not agree with the predicted ones, but investigation of this 
effect is beyond the scope of the current paper. One inter- 
esting point is that the nonlinear LDC for the quadratic law 
dips below zero at the shortest wavelengths, an effect which 
is not present in any of the LDC predictions used in this 
work. 

When combining the parameter values obtained for each 
of the ten Knutson light curves, their scatter was much larger 
than the derived uncertainties. The scatter is again rela- 
tively the largest for k (11.3a for the solutions with linear 
LD and 5.6 a for those with quadratic LD) so for the final 
result I adopt a weighted mean with uncertainties increased 
to reflect this greater scatter. Note that this effect is ex- 
pected because the individual light curves are not indepen- 
dent (they come from the same set of observations). This 
does not necessarily imply that there are any problems with 
the observations or the analysis for the quadratic-law case, 
but the much worse agreement for linear LD does indicate 
that this function is a poor representation of reality. 

I have investigated how the accuracy of the radius mea- 
surements varies with wavelength, and find (after compen- 
sating for the variation in the precision of the data) that 
the precision at the reddest wavelengths is roughly twice 
as good as at the bluest. The quality of the LDC deter- 
minations does not change much, but their importance to 
ta and rb decreases with increasing wavelength. This is in 
agreement with the general belief that redder passbands give 
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Figure 19. Variation of the linear LD law coefficient with wave- 
length of observation for the JK TEBOP solution s of the ten HST 
LCs of HD 209458 presented bv lKnutson et alj <2007ch . The un- 
certainty on each point comes from Monte Carlo simulations. The 
filled circles connected by a solid li ne show the p redicted LDCs 
for the SDSS ugriz passbands from IClareti (2004(1 . In each case 
the central wavelength of the passband is plotted; detailed con- 
clusions from these figures must take into account the passband 
response function as well as central wavelength. 
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Figure 20. Variation with wavelength of the linear and nonlin- 
ear coefficients of th e quadratic LD law f or the ten HST light 
curves of HD 209458 jKnutson et alj|2007t]) . The uncertainty on 
each point comes from Monte Carlo simulations. The filled circles 
connected by a solid l ine show the p redicted LDCs for the SDSS 
ugriz passbands from Claret (2004). The central wavelengths of 
the passbands are plotted; detailed conclusions from these figures 
must take into account the passband response function as well as 
central wavelength. 
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more precise results because the effect of LD is less impor- 
tant. 



Table 15. Stellar velocity amplitudes for the systems studied 
here. 



3. U.2 Final results for HD 209458 

The results of the three light curve analyses for HD 209458 
are given in Table[l4]and are not in good agreement. There 
is again a difference between the results which is larger than 
expected given the derived errors. This is worst for the ra- 
tio of the radii (3.3 a) but surprisingly does not affect ta or 
Tb (0.5 er). The fact that the discrepancy remains the worst 
for k is disquieting because this should be very well deter- 
mined from the observed eclipse depths. Our understanding 
of HD 209458 is good, like HD 189733, because sufficient in- 
dependent data exist for us to make a useful estimate of the 
importance of systematic errors. 

For the final result I have again adopted a weighted 
mean of the parameters from the three separate light curves 
but increased the uncertainties to account for the discrep- 
ancy between the measurements. The final parameters found 
here are in reasonable agreement with most literature val- 
ues. The orbital inclinations found in some published studies 
are not in good agreement within the errors; this can be at- 
tributed to those studies fixing LDCs at theoretical values 
and so ending up with biased best-fit parameters and un- 
derestimated errorbars. 



4 SURFACE GRAVITIES OF EACH SYSTEM 



ISouthworth et al l |2004d ) showed that it is possible to 
measure the surface gravity of the secondary component 
of an eclipsing system from only the results of a light 
curve analysis and the measured velocity amplitude of the 
primary component^, A simplified formula was given by 
ISouthworth et alj i|2007bl ): 



5b 



2vr 

-^orb 



(1 



(7) 



where K\ is the stellar velocity amplitude, and applied to 
the fourteen TEPs then known. The importance of this 
is that whilst <?b can be found directly for a transiting 
planet, measuring its other physical properties (mass, ra- 
dius, density) requires additional constraints such as are 
usually obtained f rom st ellar theory, so are not empirical. 
ISouthworth et al l (|2007bl ) found a significant correlation be- 
tween t/b and P rb, although it is not yet clear whether this 
is a real effect or arises from selection biases. A plausible 
explanation is that short-period planets are close to their 
parent star and much of their atmosphere has been evapo- 
rated off by UV and X-ray irradiation. High-energy irradi- 
ation is believed to cause close-in extra soiar planets to lose 
their atmospheres l|Lammer et alj|2003l ). 

The (?b values for the TEPs studied in this work have 
been determined from the light curve analyses and pub- 
lished measurements of Ka (Table[TSJ). These have been aug- 
mented by literature values for the other known TEPs, and 



7 This conc ept was also indepen dently discovered bv lWinn et alj 
d2007dl ') and lBeattv et alj (|2007| > - note that this d oes not neces- 
sarily require the eclipses to be total as claimed by IBeattv et alj 
|2007t) . 



System 


K A ( 


m 


s- 1 ) 


TrES-1 


115.2 


± 


6.2 


TrES-2 


181.3 


± 


2.6 


XO-1 


116.0 


± 


9.0 


WASP-1 


114 


± 


13 


WASP-2 


155 


± 


7 


T T A m T~fc 1 

HAT-P-1 


60.3 


± 


2.1 


OGLE-TR-10 


81 


± 


17 


OGLE-TR-10 


80 


± 


25 


OGLE-TR-10 


80.3 


± 


14.1 


OGLE-TR-56 


212 


± 


22 


OGLE-TR-111 


78 


± 


11 


OGLE-TR-132 


141 


± 


42 


GJ 436 


18.34 


± 


0.52 


HD 149026 


43.3 


± 


1.2 


HD 189733 


205 


± 


6 


HD 209458 


85.1 


± 


1.0 



Reference 



Alonso et al. (2004) 
O'Donovan et al. (2006) 
McCullough et al. (20061 
Collier Camerone^aL (2 
Collier Cameron ct al. (2 
Bakos et al. ('2007a') 
Bouchv et al. (2005a) 
Konacki et al. (2005) 
(adopted in this work) 
Bouchv et al. (2005a) 
Pont et al. (2004) 
Bouchv et al. (2004) 
Maness et al. (2007) 



Sato et al. (2005) 
Bouchv et al. (2005b) 
Naef et al. (2004) 



are given in Table[l6]and plotted in Fig. 1211 The correlation 
is still clearly visible but is quite weak. The error bars for 
some of the planets could be decreased by performing light 
curve analyses similar to those presented in this work. The 
fractional radii from the light curve analyses are sensitive to 
the amount of stellar irradiation intercepted by the transit- 
ing planets, but plots of ta and rb show no correlation with 

Porb • 



5 SUMMARY AND CONCLUSIONS 

At this time over thirty transiting extrasoiar planetary sys- 
tems are known, but their properties have been obtained 
using a variety of different methods. I have presented de- 
tailed light curve analyses of fourteen TEPs for which good 
light curves are available, paying particular attention to the 
importance of limb darkening (LD) and the assessment of 
random and systematic errors. In a second paper the re- 
sults of these analyses will be used to determine the physical 
properties of each system, in an homogeneous way, using the 
predictions of theoretical stellar models and concentrating 
on obtaining the most reliable values and uncertainties. 

The parameters found in my light curve analyses are 
generally in good agreement with published studies. How- 
ever, the uncertainties found in this study are often signif- 
icantly larger (by factors of up to 3). This stems from in- 
cluding several sources of error (Monte Carlo simulations for 
random errors, a residual permutation algorithm to assess 
systematic errors, and the effect of using different LD laws) , 
some of which are not incorporated into previous works. 

Although my uncertainty estimates are generally larger 
than those from other studies, solutions of different light 
curves for individual systems still do not agree to within 
these errors. This problem is clear for the two well-studied 
systems HD 189733 and HD 209458, both of which have 
three different high-quality light curves available. The worst 
discrepancy occurs for the ratio of the radii (6.7a for 
HD 189733 and 3.7a for HD 209458), which is primarily de- 
pendent on the observed transit depth. This indicates that 
systematic errors remain in the light curves which cannot 
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Table 14. Final parameters of the fit to the lBrown et al.l l|200ll ). lKnutson et al.l [|2007d ) and lRowe et al.l (|2006h light curves of HD 209458 
from the JKTEBOP analysis presented here. These have been compared to the properties of HD 209458 found in the literature. The 
uncertainties quoted for the final results from this work have been increased to yield X r e<j = 1-0- Quantities without quoted uncertainties 
were not given directly but have been calculated from other parameters which were. 
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Figure 21. Comparison between the surface gravities and orbital periods of the known transiting exoplanets. Filled circles denote the 
values obtained in this work and and open circles the systems for which values have been taken from the literature. The error bars for 
HD 209458 (P rb = 3.52 d) are smaller than the plotted symbol. A close-up of the main population of known TEPs is shown in the inset 
panel. 



be detected through analysis of only one dataset; definitive 
results for an individual TEP can only be obtained when at 
least three independent light curves are obtained and anal- 
ysed. 

I was not able to reproduce published solutions of the 
light curves for HD 149026, which contains a high-density 
planet. The three light curves of this system give quite dif- 
ferent results and it is clear that further observations are 



needed. My solution has a higher orbital inclination (corre- 
sponding to a central transit) and leads to smaller radii for 
both the star and planet, which suggests that HD 149026 b 
is even more dense than currently thought. 

An improved light curve is also needed for OGLE-TR- 
132 because the current best dataset does not cover the full 
transit. Further observations would also be useful for TrES- 
2, the only system where light curve solutions with one or 
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Table 16. Measured planetary surface gravities for the systems 
studied here (upper part of table) and for the other known TEPs 
(lower part of table) . References are entirely parenthesised if the 
surface gravity has been calculated from the planetary properties 
(Mb and iJt>) given in that work. 



system 


9b ( m s ) 


txctcrcncc 


TrES-l 


15.6±1.2 


This work 


TrES-2 


19.9±1.2 


This work 


AU-1 


15.8 ± 1.5 


This work 


WASP-l 


9.9t?; 7 


This work 


WASP-2 


19.7±2.4 


This work 


HAT-P-1 


9.05 ±0.66 


This work 


OGLE-TR-10 


5.7±1.2 


This work 


OGLE-TR-56 


22.2±6.6 


This work 


OGLE-TR-111 


11.5±2.5 


This work 


OGLE-TR-132 


15.6±6.1 


This work 


GJ 436 


12.8±1.2 


This work 


HD 149026 


„ 7+4.0 
"''-6.2 


This work 


HD 189733 


22.0 ±1.4 


This work 


HD 209458 


9.08 ±0.17 


This work 



TrES-3 

TrES-4 

XO-2 

XO-3 

WASP-3 

WASP-4 

WASP-5 

OGLE-TR-113 

OGLE-TR-182 

OGLE-TR-211 

HAT-P-2 

HAT-P-3 

HAT-P-4 

HAT-P-5 

HAT-P-6 

HD 17156 



28.4±4.2 
7.43 ±0.65 
14.9 ±1.7 
89 ±11 
23.4 ±5.9 
13.91" 
30.5±|;f 
27.6 ±4.1 

13.7H:? 

236.ltf 2 2 5 
20.4 ±3.1 

10.47 ±0.48 
16.5±1.9 
14.8 ±1.0 
11.6±5.7 



(O'Donovan et al. 2006) 
O'Donovan et al. (2007) 
(Burke et al. 20071 
fJohns-Krull et al. 2007) 
Pollacco et al. (2007) 
Wilson et al. (2008) 
Anderson et al. (2008) 
(Gillon et al. 2006) 
(Pont et al. 2007a) 
(Udalski et al. 2007) 
Loeillet et al. (2007) 
Torres et al. (2007) 
Kovacs et al. (2007) 
Bakos et al. (2007b) 
Noves et al. (2008) 
(Barbieri et al. 2007) 



two fitted LD coefficients (LDCs) do not agree well. Addi- 
tional investigations of OGLE-TR-10 would also be useful, 
because published studies of this object give widely varying 
results. 

The correlation between orb ital period and planetary 
surface gravity (first noticed by ISouthworth et all [2007b) 
remains visible when the new TEPs are included, although 
there are a few outliers with much higher gravities. Sev- 
eral of the surface gravity measurements are quite imprecise, 
with the main contribution to their uncertainty arising from 
the stellar velocity amplitude rather than the light curve 
analyses. New spectroscopic radial velocity studies for XO- 
1, WASP-l, WASP-2 and the OGLE transiting planets are 
needed to improve our understanding of the physical prop- 
erties of these systems. 



5.1 Treatment of limb darkening 

One potentially thorny problem in the study of the light 
curves of TEPs is how much effect the use of different LD 
laws has on the light curve solutions. The linear LD law is 
clearly inadequate when studying the high-quality HST light 
curves of HD 209458 and can be ruled out at the 99.99% 



level. This is in agree ment with results from th e study of 
eclipsing binary stars ijSouthworth et all l2007al) and with 
the predictions of stellar model atmospheres I Van Hammd 
119951 ; IClaretll2000l ). However, the linear law is adequate for 
most of the datasets studied in this work (particularly those 
from longer wavelengths). 

Solutions were also done using the two-coefficient 
quadratic, square-root, logarithmic and cubic LD laws. It is 
usually not possible to measure both LDCs for one star be- 
cause they are strongly correlated, but this means that one 
LDC can be fixed at a reasonable value without making a 
significant difference to the light curve solution. When both 
LDCs are included as fitted parameters, the cubic law gives 
the best-determined LDCs and the square-root law suffers 
from the biggest correlations. The quadratic and cubic laws 
generally perform well, but the usefulness of the cubic law is 
restricted as there are no theoretical predictions available. 
In general, though, there is little to choose between solu- 
tions using the four different LD laws: the solutions of all 
but the highest-quality transit light curves are not adversely 
affected by the choice of LD law. This is good news for stud- 
ies of transiting planets (but bad news for studies of stellar 
limb darkening). 

Additional solutions with all LDCs fixed were calculated 
for TrES-l, TrES-2 and XO-1. In general these give results 
in (mostly) acceptable agreement with other solutions, but 
with substantially smaller errorbars. It is therefore clear that 
fixing coefficients at theoretical values not only allows biases 
to creep in to the results, but also makes the uncertainties 
too optimistic. 

Despite the previous point, a comparison between mea- 
sured LDCs and those predicted using stellar model atmo- 
spheres is a promising way of improving our knowledge of 
stellar atmospheres. It is known that current stellar model 
atmospheres do not reproduce the LD of stars very well 
(A. Claret, private communication). Whilst a detailed study 
is beyond the scope of the present work, a limited com- 
parison was carried out between the LDCs measured from 
light curves of HD 209458 in ten passbands spread through 
the optical wavelength range. Predicted LDCs were unable 
to match the measured coefficien ts for either the line ar or 
quadratic laws. A similar study bv lTinglev et al. I (120061) also 
found a disagreement between the observed and predicted 
colour dependence of limb-darkened transit light curves. 



5.2 Looking ahead 

Much of our understanding of extrasolar planets rests on the 
transiting systems, and this work confirms that these can 
be studied to high accuracy. New light curves of HD 149026, 
TrES-2, OGLE-TR-10, OGLE-TR-56, OGLE-TR-132 and 
GJ436 would help to improve our understanding further, 
as would additional radial velocity studies of XO-1, WASP- 
1, WASP-2 and the OGLE transiting planets. In Paper II 
(in preparation) the light curve analysis results from the 
current paper will be used to determine homogeneous sets of 
physical properties for the fourteen TEPs studied here, with 
the aim of improving the prospects for statistical studies of 
these and other systems. 
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APPENDIX A: RESULTS OF THE LIGHT 
CURVE ANALYSES 

The tables in this section contain the detailed results of the 
modelling of each light curve. Note that whilst all the re- 
sults are best fits to the relevant data, some parameters are 
unphysical (for example the limb darkening coefficients im- 
ply that the limb of the star produces a negative amount of 
light). In these cases the unphysical results have not been 
used but are retained in the table for completeness. 
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of LD. The parameter symbols are as given in Section[2]2] The transit midpoint To is expressed as HJD — 2 453 000. The perturbations 
of the LD coefficients are discussed in Section fe.SI 



Linear LD law Quadratic LD law Square-root LD law Logarithmic LD law Cubic LD law 

All LD coefficients fixed 

, 1082+ - 0018 1078+ ' 0015 1078+ ' 0016 1077+ ' 0016 

?A + 'b u - luoz -0.0004 U ' iU ' 8 -0.0004 U - iU ' 8 -0.0005 U - iU ' '-0.0004 

k - 13444 to.ooo38 0.13568 ±0.00039 0.13533 ± 0.00040 0.13566 ± 0.00041 

• (dag.) 89.85±g;J* 89.96t°;°* 89.91±°;°? 89.99^°;^ 

M A 0.53 fixed 0.28 fixed 0.14 fixed 0.61 fixed 

v A 0.29 fixed 0.52 fixed 0.26 fixed 

T 898.87336 ±0.00009 898.87336 ± 0.00008 898.87336 ±0.00008 898.87336 ± 0.00008 



r A 0-0954I—S U.U949I-— I U - U95U Io.5oi3 U - Uy4 8_a5oi4 

n ni98q+ 00007 n ni98R+ 00007 n m ooc+o.oooos n nl 9 o fi +o.oooos 

r b 0.01283_ 00027 0.01288_ 0Q020 0.01285_ 00023 0.01286__ 0Q023 

a (mmag) 1.6785 1.6420 1.6478 1.6425 

X 2 rod 1.0544 1.0085 1.0157 1.0091 



Fitting for one LD coefficient 

r A + r b 0.1097±0.0015 0.1087 ± 0.0016 0.1094 ± 0.0019 0.1091 ± 0.0017 0.1091 ± 0.0017 

k 0.13802 ±0.00073 0.13654 ± 0.00070 0.13729 ± 0.00073 0.13693 ± 0.00075 0.13691 ± 0.00070 

i (deg.) 88.69±0.66 89.18±0.65 88.86±0.54 88.99±0.65 88.99±0.64 

u A 0.368 ±0.029 0.245 ±0.027 0.071 ±0.030 0.558 ±0.037 0.331 ±0.024 

v A 0.29 perturbed 0.52 perturbed 0.26 perturbed 0.20 perturbed 

T 898.87336 ±0.00010 898.87336 ± 0.00009 898.87336 ± 0.00008 898.87336 ± 0.00008 898.87336 ± 0.00008 



r A 0.0951 ±0.0011 0.0956 ± 0.0013 0.0962 ± 0.0016 0.0959 ± 0.0014 0.0960 ± 0.0014 

r b 0.01293 ±0.00020 0.01306 ± 0.00024 0.01321 ± 0.00029 0.01314 ± 0.00026 0.01314 ± 0.00025 

a (mmag) 1.6386 1.6408 1.6392 1.6397 1.6398 

X 2 rod 1.0052 1.0078 1.0059 1.0065 1.0066 



Fitting for two LD coefficients 

r A + r b 0.1097 ±0.0020 0.1097 ± 0.0018 0.1098 ± 0.0019 0.1096 ± 0.0018 

k 0.1380±0.0010 0.1382±0.0012 0.1382 ± 0.0013 0.1382 ± 0.0012 

i (deg.) 88.70±0.65 88.69±0.52 88.66±0.53 88.67±0.48 

u A 0.37 ±0.08 0.45 ±0.41 0.34 ±0.18 0.38 ±0.05 

v A 0.00 ±0.17 -0.13 ±0.71 -0.04 ±0.24 -0.05 ±0.19 
T 898.87336 ±0.00008 898.87336 ± 0.00008 898.87336 ± 0.00008 898.87336 ± 0.00009 



r A 0.0964 ±0.0016 0.0964 ± 0.0015 0.0964 ± 0.0016 0.0963 ± 0.0015 

r b 0.01330 ±0.00031 0.01332 ± 0.00032 0.01332 ± 0.00032 0.01332 ± 0.00030 

o- (mmag) 1.6387 1.6386 1.6386 1.6386 

X 2 rcd 1.0061 1.0060 1.0060 1.0060 
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Table A2. Parameters of the JKTEBOP best fits of the TrES-2 light curve (z band, N ohs = 1033) from lHolman et aD l|2007bl 'l for different 
treatments of LD. The parameter symbols are as given in Section ^. 21 The transit midpoint To is expressed as HJD — 2 453 000. For each 
part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


All LD coefficients fixed 


?A + r b 


0.1468 ±0.0020 


0.1475 ±0.0020 


0.1479 ±0.0020 


0.1478 ±0.0021 




k 


0.12681 ±0.00086 


0.12537 ±0.00086 


0.12537 ±0.00087 


0.12532 ±0.00089 




i (deg.) 


83.61±0.12 


83.58±0.12 


83.54 ±0.12 


83.55 ±0.12 




"A 


0.49 fixed 


0.22 fixed 


0.06 fixed 


0.57 fixed 








0.31 fixed 


0.56 fixed 


0.27 fixed 




± 


989 75282 ± 0001 7 


989 75298 + 0001 8 


989 75279 + 0001 7 


989 75279 + 0001 8 




TA 


0.1303 ±0.0017 


0.1310 ±0.0017 


0.1314 ±0.0018 


0.1313 ±0.0017 






0.01653 ±0.00028 


0.01643 ±0.00030 


0.01647 ±0.00030 


0.01645 ±0.00031 




a (mmag) 


1.4941 


1.4940 


1.4942 


1.4941 




V 2 

* red 


1.0050 


1.0048 


1.0052 


1.0050 




Fitting for 


one LD coefficient 










TA + r b 


0.1453 ±0.0035 


0.1462 ±0.0036 


0.1461 ±0.0038 


0.1462 ±0.0037 


0.1455 ±0.0037 


k 


0.1281 ±0.0025 


0.1263 ±0.0028 


0.1269 ±0.0026 


0.1265 ±0.0029 


0.1274 ±0.0026 


i (deg.) 


83.75 ±0.31 


83.69 ±0.37 


83.71 ±0.37 


83.70 ±0.38 


83.75 ±0.36 


«A 


0.61 ±0.22 


0.33 ±0.30 


0.21 ±0.27 


0.69 ±0.32 


0.55 ±0.26 


va 




0.31 perturbed 


0.56 perturbed 


0.27 perturbed 


0.10 perturbed 


To 


989.75279 ±0.00018 


989.75279 ±0.00018 


989.75279 ±0.00017 


989.75279 ±0.00018 


989.75279 ±0.00018 


TA 


0.1288 ±0.0032 


0.1298 ±0.0033 


0.1296 ±0.0034 


0.1298 ±0.0034 


0.1291 ±0.0033 


r b 


0.01650 ±0.00040 


0.01640 ±0.00038 


0.01645 ±0.00041 


0.01643 ±0.00041 


0.01645 ±0.00039 


a (mmag) 


1.4939 


1.4939 


1.4939 


1.4939 


1.4939 


Y 2 

A red 


1.0057 


1.0057 


1.0057 


1.0057 


1.0057 


Fitting for two LD coefficients 


TA + r b 




0.1431 ±0.0053 


0.1451 ±0.0049 


0.1428 ±0.0064 


0.1370 ±0.0059 


k 




0.1316 ±0.0073 


0.1282 ±0.0048 


0.1322 ±0.0086 


0.1409 ±0.0090 


i (deg.) 




83.91 ±0.41 


83.77 ±0.39 


83.91 ±0.44 


84.29 ±0.39 


«A 




1.07 ±0.77 


0.62 ±1.27 


0.41 ±0.58 


1.46 ±0.44 


VA 




-0.50 ±0.95 


-0.01 ±1.83 


-0.60 ±1.21 


-1.17±0.68 


To 




989.75280 ±0.0018 


989.75279 ±0.00016 


989.75279 ±0.00018 


989.75278 ±0.00018 


TA 




0.1264 ±0.0050 


0.1286 ±0.0047 


0.1261 ±0.0065 


0.1201 ±0.0059 


r b 




0.01664 ±0.00050 


0.01649 ±0.00043 


0.01667 ±0.00052 


0.01692 ±0.00047 


a (mmag) 




1.4938 


1.4939 


1.4937 


1.4937 


Y 2 

red 




1.0065 


1.0066 


1.0064 


1.0061 
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Table A3. Parameters of the JKTEBOP best fits of the XO-1 light curve from the FLWO 1.2m telescope l|Holman et aljlioolih (z band, 
N obs = 821) for different treatment of LD. The parameter symbols are as given in Section ^, 21 The transit midpoint To is expressed as 
HJD — 2 453 000. For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


All LD coefficients fixed 


?A + r b 


0.1003 ±0.0019 


0.1004 ±0.0017 


0.1012 ±0.0019 


0.1010 ±0.0018 




k 


0.1317 ±0.0007 


0.1327 ±0.0005 


0.1331 ±0.0006 


0.1330 ±0.0006 




i (deg.) 


88.94 ±0.57 


88.87 ±0.42 


88.66 ±0.37 


88.71 ±0.38 




"A 


0.49 fixed 


0.23 fixed 


0.07 fixed 


0.57 fixed 








0.31 fixed 


0.55 fixed 


0.27 fixed 




± 


875 92236 ± 0001 


875 92236 + 0001 1 


875 92235 + 0001 1 


875 92236 + 0001 1 




TA 


0.0886 ±0.0016 


0.0886 ±0.0015 


0.0893 ±0.0016 


0.0892 ±0.0015 






0.01167 ±0.00027 


0.01176 ±0.00023 


0.01189 ±0.00026 


0.01186 ±0.00025 




a (mmag) 


1.5081 


1.4816 


1.4841 


1.4826 




V 2 

* red 


1.0568 


1.0179 


1.0217 


1.0195 




Fitting for 


one LD coefficient 










TA + r b 


0.1029 ±0.0017 


0.1004 ±0.0019 


0.1022 ±0.0018 


0.1017 ±0.0020 


0.1025 ±0.0018 


k 


0.1347 ±0.0007 


0.1327 ±0.0009 


0.1338 ±0.0008 


0.1334 ±0.0009 


0.1342 ±0.0008 


i (deg.) 


88.26 ±0.28 


88.87 ±0.54 


88.44 ±0.33 


88.55 ±0.38 


88.36 ±0.30 


«A 


0.358 ±0.025 


0.230 ±0.031 


0.042 ±0.030 


0.557 ±0.036 


0.430 ±0.033 


va 




0.31 perturbed 


0.55 perturbed 


0.27 perturbed 


0.10 perturbed 


To 


883.80543 ±0.00011 


883.80542 ±0.00010 


883.80542 ±0.00010 


883.80542 ±0.00011 


883.80542 ± 0.00010 


TA 


0.0906 ±0.0015 


0.0886 ±0.0017 


0.0901 ±0.0016 


0.0898 ±0.0017 


0.0903 ±0.0015 


r b 


0.01221 ±0.00025 


0.01176 ±0.00029 


0.01206 ±0.00027 


0.01198 ±0.00029 


0.01213 ±0.00026 


a (mmag) 


1.4849 


1.4832 


1.4824 


1.4824 


1.4838 


Y 2 

A red 


1.0244 


1.0218 


1.0206 


1.0206 


1.0228 


Fitting for two LD coefficients 


TA + r b 




0.1001 ±0.0019 


0.1007 ±0.0020 


0.1004 ±0.0019 


0.1003 ±0.0019 


k 




0.1320±0.0013 


0.1310 ±0.0021 


0.1316 ±0.0017 


0.1314 ±0.0016 


i (deg.) 




89.03 ±0.72 


89.06 ±0.74 


89.05 ±0.73 


89.08 ±0.72 


«A 




0.17±0.08 


-0.86 ±0.63 


0.84 ±0.25 


0.86 ±0.23 


VA 




0.46 ±0.19 


2.17±1.13 


0.63 ±0.31 


0.65 ±0.30 


To 




883.80542 ±0.00011 


883.80542 ±0.00011 


883.80542 ±0.00010 


883.80542 ±0.00011 


TA 




0.0885 ±0.0016 


0.0890 ±0.0018 


0.0887 ±0.0017 


0.0887 ±0.0016 


r b 




0.01167 ±0.00031 


0.01166 ±0.00031 


0.01167 ±0.00030 


0.01165 ±0.00029 


a (mmag) 




1.4809 


1.4807 


1.4808 


1.4808 


Y 2 

red 




1.0195 


1.0192 


1.0193 


1.0193 
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Table A4. Parameters of the JKTEBOP best fits of the WASP-1 light curve from lShporer et al ] (|2007t l (I band, N ohs = 583) for different 
treatment of LD. The parameter symbols are as given in Section[2]2] The transit midpoint To is expressed as HJD — 2 454 000. For each 
part of the table the upper parameters are fitted quantities and the lower parameters arc derived quantities. 







Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










r\ + r h 




0.2036 ±0.0104 


0.1986 ±0.0102 


0.2032 ±0.0106 


0.1978 ±0.0096 


0.2027 ±0.0114 


k 




0.1042 ±0.0011 


0.1031 ±0.0012 


0.1038 ±0.0012 


0.1025 ±0.0013 


0.1037 ±0.0012 


i (deg.) 




86.0±1.6 


87.0 ±2.4 


86.1 ±2.2 


87.3 ±2.3 


86.2±2.3 


u A 




0.215 ±0.070 


0.043 ±0.082 


-0.114 ±0.079 


0.644 ±0.088 


0.158 ±0.074 


VA 






0.35 perturbed 


0.55 perturbed 


0.60 perturbed 


0.20 perturbed 


T 




13.3127 ±0.0004 


13.3127 ±0.0004 


13.3127 ±0.0004 


13.3127 ±0.0004 


13.3127 ±0.0004 


TA 




0.1844 ±0.0093 


0.1800 ±0.0092 


0.1841 ±0.0095 


0.1794 ±0.0086 


0.1837 ±0.0102 


rb 




0.0192 ±0.0011 


0.0186 ±0.0011 


0.0191 ±0.0011 


0.0184 ±0.0010 


0.0190 ±0.0012 


a (mmaj 





2.1653 


2.1664 


2.1655 


2.1665 


2.1656 


Y 2 

^ red 




1.0016 


1.0032 


1.0020 


1.0034 


1.0022 


Fitting for two LD coefficients 


rA + r h 






0.2038 ±0.0106 


0.2035 ±0.0106 


0.2035 ±0.0109 


0.2034 ±0.0110 


k 






0.1044 ±0.0017 


0.1044 ±0.0019 


0.1044 ±0.0019 


0.1043 ±0.0019 


i (deg.) 






85.9±1.9 


86.0 ±1.2 


86.0±1.9 


86.0 ±2.2 


u A 






0.29 ±0.32 


0.41 ±1.17 


0.14 ±0.44 


0.23 ±0.17 


VA 






-0.13 ±0.54 


-0.32 ±1.90 


-0.12 ±0.68 


-0.05 ±0.54 


To 






13.3127 ±0.0004 


13.3127 ±0.0004 


13.3127 ±0.0004 


13.3127 ±0.0004 


?"A 






0.1845 ±0.0094 


0.1842 ±0.0094 


0.1843 ±0.0098 


0.1842 ±0.0099 


rb 






0.0193 ±0.0013 


0.0192 ±0.0013 


0.0192 ±0.0013 


0.0192 ±0.0013 


cr (mmaj 







2.1653 


2.1653 


2.1653 


2.1653 


Y 2 

A red 






1.0032 


1.0033 


1.0033 


1.0033 



Table A5. Parameters of the JKTEBOP best fits of the WASP-1 light curve from ICharbonneau et al.l (120071 ') (z band, iV obs = 657) for 
different treatment of LD. The parameter symbols are as given in Section ^, 21 The transit midpoint To is expressed as HJD — 2 450 000. 
For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










r A + T'b 


0.1941 ±0.0065 


0.1918 ±0.0052 


0.1922 ±0.0057 


0.1912 ±0.0047 


0.1919 ±0.0052 


k 


0.10479 ±0.00097 


0.10405 ±0.00093 


0.10429 ±0.00093 


0.10411 ±0.00092 


0.10402 ±0.00089 


i (deg.) 


88.0±1.9 


89.3±1.6 


88.8 ±1.7 


89.6±1.6 


89.2±1.6 


UA 


0.211 ±0.54 


0.055 ±0.060 


-0.141 ±0.058 


0.408 ±0.059 


0.143 ±0.061 


VA 




0.32 perturbed 


0.60 perturbed 


0.28 perturbed 


0.28 perturbed 


To 


5.7520 ±0.0003 


5.7520 ±0.0003 


5.7520 ±0.0003 


5.7520 ±0.0003 


5.7520 ±0.0003 


VA 


0.1757 ±0.0058 


0.1737 ±0.0047 


0.1740 ±0.0050 


0.1732 ±0.0042 


0.1738 ±0.0046 


rb 


0.01841 ±0.00071 


0.01807 ±0.00056 


0.01815 ±0.00062 


0.01803 ±0.00050 


0.01808 ±0.00054 


a (mmag) 


1.8084 


1.8096 


1.8086 


1.8089 


1.8088 


Y 2 

A red 


0.9616 


0.9629 


0.9618 


0.9621 


0.9621 


Fitting for two LD coefficients 


+ r h 




0.1922 ±0.0057 


0.1916 ±0.0057 


0.1921 ±0.0056 


0.1915 ±0.0055 


k 




0.10465 ±0.00114 


0.10447 ±0.00131 


0.10456 ±0.00128 


0.10446 ±0.00134 


i (deg.) 




88.6±1.8 


88.9 ±1.7 


88.7±1.8 


88.9±1.7 


u A 




0.21 ±0.20 


0.06 ±0.88 


0.25 ±0.33 


0.19±0.12 


VA 




0.01 ±0.36 


0.27 ±1.47 


0.06 ±0.47 


0.07 ±0.41 


To 




5.7520 ±0.0003 


5.7520 ±0.0003 


5.7520 ±0.0003 


5.7520 ±0.0003 


VA 




0.1740 ±0.0051 


0.1735 ±0.0051 


0.0050 ±0.0050 


0.1734 ±0.0050 


rb 




0.01820 ±0.00064 


0.01812 ±0.00065 


0.00062 ±0.00062 


0.01811 ±0.00059 


a (mmag) 




1.8084 


1.8085 


1.8084 


1.8084 


^ red 




0.9631 


0.9632 


0.9631 


0.9631 
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Table A6. Parameters of the JKTEBOP best fits of the WASP-2 light curve from [ctiarbonne au et alj ((2007) (z band, N obs = 426) for 
different treatment of LD. The parameter symbols are as given in Section[2]2] The transit midpoint To is expressed as HJD — 2 454 000. 
For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










fA + r b 


0.1403 ±0.0058 


0.1409 ±0.0059 


0.1408 ±0.0060 


0.1409 ±0.0060 


0.1409 ±0.0060 


k 


0.1322 ±0.0014 


0.1310 ±0.0015 


0.1315 ±0.0015 


0.1313 ±0.0015 


0.1313 ±0.0016 


i (deg.) 


84.86 ±0.45 


84.83 ±0.46 


84.84 ±0.47 


84.83 ±0.48 


84.83 ±0.47 


u A 


0.476 ±0.147 


0.240 ±0.171 


0.128 ±0.163 


0.594 ±0.169 


0.369 ±0.167 


VA 




0.30 perturbed 


0.53 perturbed 


0.25 perturbed 


0.20 perturbed 


T 


8.73204 ±0.00021 


8.73205 ±0.00022 


8.73205 ±0.00022 


8.73205 ±0.00022 


8.73205 ±0.00021 


VA 


0.1239 ±0.0050 


0.1246 ±0.0050 


0.1244 ±0.0051 


0.1245 ±0.0052 


0.1246 ±0.0052 


rb 


0.01637 ±0.00079 


0.01632 ±0.00079 


0.01636 ±0.00083 


0.01635 ±0.00081 


0.01635 ±0.00081 


a (mmag) 


1.8891 


1.8891 


1.8891 


1.8891 


1.8891 


Y 2 

A red 


0.9988 


0.9989 


0.9989 


0.9989 


0.9989 


Fitting for two LD coefficients 


rA + ry 




0.1405 ±0.0074 


0.1404 ±0.0066 


0.1404 ±0.0071 


0.1404 ±0.0075 


k 




0.1318 ±0.0099 


0.1321 ±0.0082 


0.1322 ±0.0093 


0.1321 ±0.0090 


i (deg.) 




84.84 ±0.47 


84.85 ±0.41 


84.85 ±0.41 


84.85 ±0.46 


MA 




0.40 ±1.63 


0.43 ±3.24 


0.47±1.19 


0.46 ±0.67 


VA 




0.10 ±2.20 


0.06 ±5.11 


0.00 ±2.22 


0.02 ±1.54 


To 




8.73205 ±0.00021 


8.73205 ±8.73205 


8.73204 ±0.00021 


8.73204 ±0.00021 


VA 




0.1241 ±0.0071 


0.1240 ±0.0063 


0.1240 ±0.0069 


0.1240 ±0.0072 


rb 




0.01636 ±0.00097 


0.01639 ±0.00091 


0.01639 ±0.00088 


0.01638 ±0.00095 


a (mmag) 




1.8891 


1.8891 


1.8891 


1.8891 


Y 2 

A rod 




1.0012 


1.0012 


1.0012 


1.0012 



Table A7. Parameters of the JKTEBOP best fits of the FLWO light curve of HAT-P-1 from lWinn et al.1 d2007d l for different treatments 
of LD. The parameter symbols are as given in Section ^, 2 1 The transit midpoint To is expressed as HJD — 2 453 000. For each part of the 
table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










^A + r h 


0.1025 ±0.0029 


0.1026 ±0.0031 


0.1026 ±0.0030 


0.1026 ±0.0031 


0.1026 ±0.0029 


k 


0.11328 ±0.00082 


0.11217 ±0.00092 


0.1126 ±0.00092 


0.1124 ±0.00087 


0.1128 ±0.00091 


i (deg.) 


86.36 ±0.22 


86.38 ±0.26 


86.37 ±0.25 


86.37±0.25 


86.37 ±0.23 


MA 


0.485 ± 0.063 


0.273 ±0.076 


0.117 ±0.070 


0.629 ±0.075 


0.439 ±0.067 


^A 




0.30 perturbed 


0.58 perturbed 


0.28 perturbed 


0.10 perturbed 


To 


997.79205 ±0.00030 


997.79205 ±0.00030 


997.79205 ±0.00026 


997.79202 ±0.00028 


997.79206 ±0.00030 


rA 


0.0920 ±0.0025 


0.0922 ±0.0028 


0.0922 ±0.0026 


0.0923 ±0.0027 


0.0922 ±0.0026 




0.01043 ±0.00035 


0.01038 ±0.00038 


0.01038 ±0.00037 


0.01037 ±0.00037 


0.01040 ±0.00036 


a (mmag) 


2.1309 


2.1320 


2.1316 


2.1318 


2.1314 


Fitting for two LD coefficients 


rA + r b 




0.1024 ±0.0018 


0.1027 ±0.0017 


0.1026 ±0.0018 


0.1028 ±0.0018 


k 




0.1221 ±0.0037 


0.1208 ±0.0028 


0.1212 ±0.0031 


0.1210 ±0.0031 


i (deg.) 




86.27±0.13 


86.22 ±0.12 


86.24 ±0.12 


86.22 ±0.13 


ua 




1.79 ±0.39 


3.62 ±0.89 


-0.31 ±0.21 


1.18 ±0.23 


VA 




-1.66 ±0.46 


-4.73 ±1.31 


-1.73 ±0.49 


-1.29 ±0.37 


To 




997.79206 ±0.00029 


997.79208 ±0.00028 


997.79208 ±0.00028 


997.79208 ±0.00028 


rA 




0.0913 ±0.0015 


0.0916 ±0.0015 


0.0915 ±0.0016 


0.0917 ±0.0016 


rb 




0.01115 ±0.00039 


0.01106 ±0.00031 


0.01109 ±0.00034 


0.01109 ±0.00036 


cr (mmag) 




2.1248 


2.1248 


2.1249 


2.1248 
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Table A8. Parameters of the JKTEBOP best fits of the Lick light curve of HAT-P-1 from|w inn et al for different treatments of 

LD. The parameter symbols are as given in Section[2]2] The transit midpoint To is expressed as H JD — 2 453 000. For each part of the 
table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










fA + n> 


0.1065 ±0.0028 


0.1055 ±0.0027 


0.1056 ±0.0030 


0.1056 ±0.0029 


0.1056 ±0.0030 


k 


0.11108 ±0.00062 


0.10992 ±0.00059 


0.11031 ±0.00066 


0.11012 ±0.00068 


0.11054 ±0.00066 


i (deg.) 


86.01 ±0.22 


86.11 ±0.22 


86.10 ±0.24 


86.10 ±0.22 


86.09 ±0.23 


u A 


0.339 ±0.076 


0.124 ±0.081 


-0.022 ±0.081 


0.486 ±0.083 


0.304 ±0.079 


VA 




0.30 perturbed 


0.58 perturbed 


0.28 perturbed 


0.10 perturbed 


To 


997.79357 ±0.00028 


997.79355 ±0.00026 


997.79355 ±0.00026 


997.79357 ±0.00028 


997.79356 ±0.00026 


TA 


0.0959 ±0.0024 


0.0950 ±0.0024 


0.0951 ±0.0027 


0.0952 ±0.0026 


0.0951 ±0.0027 


rb 


0.01065 ±0.00031 


0.01045 ±0.00031 


0.01049 ±0.00034 


0.01048 ±0.00032 


0.01051 ±0.00034 


a (mmag) 


1.9210 


1.9207 


1.9208 


1.9207 


1.9209 


Fitting for two LD coefficients 


TA + T b 




0.1019 ±0.0051 


0.1079 ±0.0055 


0.1068 ±0.0052 


0.1064 ±0.0054 


k 




0.1023 ±0.0044 


0.1029 ±0.0027 


0.1021 ±0.0034 


0.1024 ±0.0039 


i (deg.) 




86.68 ±0.40 


86.27 ±0.35 


86.31 ±0.32 


86.41 ±0.36 


u A 




-1.20 ±1.12 


-4.34 ±2.03 


2.25 ±0.95 


-0.26 ±0.41 


VA 




2.63 ±1.79 


7.78 ±3.34 


2.93 ±1.51 


2.13±1.16 


To 




997.79359 ±0.00026 


997.79360 ±0.00028 


997.79358 ±0.00024 


997.79360 ±0.00027 


TA 




0.0924 ±0.0047 


0.0978 ±0.0050 


0.0970 ±0.0049 


0.0965 ±0.0050 


rb 




0.00945 ±0.00051 


0.01007 ±0.00042 


0.01000 ±0.00039 


0.00988 ±0.00046 


cr (mmag) 




1.9190 


1.9194 


1.9194 


1.9190 



Table A9. Parameters of the JKTEBOP best fits of the OGLE-TR-10 light curve from iHolman et al.l |(2007a) (I band, A obs = 388) for 
different treatment of LD. The parameter symbols are as given in Section|2]2] The transit midpoint To is expressed as HJD — 2 453 00. 
For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 












0.1469 ±0.0058 


0.1468 ±0.0064 


0.1469 ±0.0060 


0.1468 ±0.0052 


0.1469 ±0.0061 


k 


0.1026 ±0.0010 


0.1021 ±0.0011 


0.1023 ±0.0011 


0.1021 ±0.0011 


0.1022 ±0.0011 


i (deg.) 


85.98 ±0.59 


86.05 ±0.68 


86.02 ±0.61 


86.04 ±0.65 


86.02 ±0.62 


UA 


0.223 ±0.063 


0.011 ±0.075 


-0.130 ±0.070 


0.396 ±0.071 


0.143 ±0.074 


VA 




0.32 perturbed 


0.60 perturbed 


0.30 perturbed 


0.20 perturbed 


P (days) 


3.101257 ±0.000005 


3.101258 ±0.000005 


3.101257 ±0.000005 


3.101257 ±0.000005 


3.101258 ±0.000005 


To 


583.64596 ±0.00028 


583.64593 ±0.00028 


583.64594 ±0.00028 


583.64593 ±0.00027 


583.64594 ±0.00029 


VA 


0.1332 ±0.0053 


0.1332 ±0.0057 


0.1333 ±0.0054 


0.1332 ±0.0055 


0.1333 ±0.0054 


rb 


0.01368 ±0.00062 


0.01359 ±0.00068 


0.01363 ±0.00066 


0.01361 ±0.00068 


0.01362 ±0.00067 


a (mmag) 


1.6660 


1.6691 


1.6675 


1.6684 


1.6683 


^ red 


0.9747 


0.9745 


0.9744 


0.9745 


0.9745 


Fitting for two LD coefficients 


ta + r b 




0.1466 ±0.0070 


0.1470 ±0.0067 


0.1465 ±0.0070 


0.1461 ±0.0067 


k 




0.1022 ±0.0019 


0.1022 ±0.0023 


0.1022 ±0.0020 


0.1022 ±0.0021 


i (deg.) 




86.04 ±0.79 


86.01 ±0.77 


86.06 ±0.80 


86.09 ±0.80 


u A 




0.07 ±0.38 


-0.20 ±1.45 


0.36 ±0.46 


0.17±0.20 


VA 




0.22 ±0.61 


0.67 ±2.37 


0.24 ±0.74 


0.15±0.56 


To 




583.64594 ±0.00028 


583.64594 ±0.00030 


583.64594 ±0.00029 


583.64595 ±0.00027 


TA 




0.1330 ±0.0062 


0.1334 ±0.0059 


0.1329 ±0.0061 


0.1326 ±0.0059 


rb 




0.01360 ±0.00082 


0.01364 ±0.00079 


0.01358 ±0.00081 


0.01355 ±0.00082 


a (mmag) 




1.6679 


1.6677 


1.6676 


1.6673 


red 




0.9768 


0.9770 


0.9769 


0.9769 



© 0000 RAS, MNRAS 000, 000-000 



32 John Southworth 



Table A10. Parameters of the JKTEBOP best fits of the VLT V-band light curve of OGLE-TR-10 (AT ob8 = 78) from lPont et alJ l|2007ch for 
different treatment of LD. The parameter symbols are as given in Section[2]2] The transit midpoint To is expressed as HJD — 2 453 000. 
For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










TA + r b 


0.1881 ±0.0080 


0.1885 ±0.0098 


0.1880 ±0.0088 


0.1882 ±0.0089 


0.1873 ±0.0092 


k 


0.1183 ±0.0014 


0.1174±0.0017 


0.1175 ±0.0015 


0.1176 ±0.0015 


0.1172 ±0.0016 


i (deg.) 


82.96 ±0.65 


82.95 ±0.79 


83.00 ±0.74 


82.98 ±0.74 


83.07 ±0.76 


"A 


0.48 ±0.14 


0.29±0.19 


0.13±0.16 


0.58 ±0.17 


0.40 ±0.17 






0.25 perturbed 


0.55 perturbed 


0.20 perturbed 


0.20 perturbed 


T 


890.67095 ±0.00047 


890.67095 ±0.00049 


890.67095 ±0.00047 


890.67095 ±0.00045 


890.67096 ±0.00047 


TA 


0.1682 ±0.0070 


0.1687 ±0.0086 


0.1682 ±0.0076 


0.1684 ±0.0078 


0.1676 ±0.0080 


rb 


0.0199 ±0.0010 


0.01980 ±0.0012 


0.0198 ±0.0011 


0.0198 ±0.0011 


0.0196 ±0.0011 


cr (mmag) 


0.9756 


0.9774 


0.9771 


0.9771 


0.9777 


Y 2 

A red 


1.8409 


1.8470 


1.8459 


1.8459 


1.8480 


Fitting for two LD coefficients 


rA + r h 




0.1835 ±0.0053 


0.1837 ±0.0049 


0.1838 ±0.0052 


0.1841 ±0.0052 


k 




0.1289 ±0.0107 


0.1274 ±0.0085 


0.1278 ±0.0096 


0.1274 ±0.0091 


i (deg.) 




83.11 ±0.40 


83.03 ±0.40 


83.05 ±0.41 


83.03 ±0.41 


MA 




1.95 ±0.90 


4.04 ±2.03 


-0.39 ±0.49 


1.27±0.54 


^A 




11.84 ±1.10 


-5.34 ±2.98 


-1.93 ±1.15 


-1.43 ±0.92 


To 




890.67098 ±0.00048 


890.67097 ±0.00044 


890.67098 ±0.00048 


890.67098 ±0.00044 


TA 




0.1626 ±0.0051 


0.1629 ±0.0044 


0.1629 ±0.0046 


0.1633 ±0.0048 


rb 




0.0210 ±0.0017 


0.0208 ±0.0014 


0.0208 ±0.0016 


0.0208 ±0.0015 


cr (mmag) 




0.9666 


0.9659 


0.9662 


0.9662 


Y 2 

A rod 




1.8336 


1.8310 


1.8321 


1.8323 



Table All. Parameters of the JKTEBOP best fits of the VLT R-band light curve of OGLE-TR-10 (7V obs = 67) from lPont et alJ ll2007ch for 
different treatment of LD. The parameter symbols are as given in Section ^, 21 The transit midpoint To is expressed as HJD — 2 453 000. 
For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










rA + r h 


0.1636 ±0.0093 


0.1634 ±0.0095 


0.1634 ±0.0096 


0.1633 ±0.0093 


0.1633 ±0.0092 


k 


0.11764 ±0.00067 


0.11720 ±0.00068 


0.11722 ±0.00069 


0.11710 ±0.00069 


0.1174 ±0.00067 


i (deg.) 


84.61 ±0.74 


84.64 ±0.76 


84.64 ±0.75 


84.65 ±0.74 


84.64 ±0.72 


MA 


0.00 ±0.11 


-0.19 ±0.13 


-0.35 ±0.12 


0.13 ±0.12 


-0.05 ±0.11 


^A 




0.25 perturbed 


0.55 perturbed 


0.25 perturbed 


0.10 perturbed 


To 


890.67095 ±0.00044 


890.67095 ± 0.00047 


890.67095 ±0.00043 


890.67095 ± 0.00043 


890.67095 ±0.00044 


TA 


0.1463 ±0.0083 


0.1463 ±0.00086 


0.1462 ±0.0085 


0.1462 ±0.0084 


0.1462 ±0.0082 


rb 


0.0172 ±0.0010 


0.0171 ±0.0010 


0.0171 ±0.0010 


0.0181 ±0.0010 


0.0172 ±0.0010 


a (mmag) 


0.7155 


0.7149 


0.7150 


0.7149 


0.7151 


^ red 


1.6819 


1.6790 


1.6797 


1.6792 


1.6802 


Fitting for two LD coefficients 


rA + r h 




0.1402 ±0.0040 


0.1575 ±0.0110 


0.1536 ±0.0077 


0.1565 ±0.0110 


k 




0.1015 ±0.0027 


0.1125 ±0.0055 


0.1050 ±0.0018 


0.1141 ±0.0059 


i (deg.) 




88.98 ±1.13 


85.54 ±1.02 


86.64 ±0.97 


85.48 ±0.99 


M A 




-0.54 ±0.22 


-3.54 ±3.53 


3.04 ±0.25 


-0.30 ±0.50 


UA 




2.56 ±0.43 


6.00 ±5.89 


3.92 ±0.48 


1.13 ±1.70 


T 




890.67093 ±0.00046 


890.67098 ±0.00048 


890.67100 ±0.00045 


890.67097 ±0.00048 


rA 




0.1272 ±0.0036 


0.1416 ±0.0099 


0.1390 ±0.0068 


0.1405 ±0.0101 


rb 




0.0129 ±0.00047 


0.0159 ±0.0013 


0.0146 ±0.00090 


0.0160 ±0.0012 


cr (mmag) 




0.6875 


0.7142 


0.7131 


0.7142 


X red 




1.5941 


1.7020 


1.6994 


1.7020 
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Table A12. Parameters of the JKTEBOP best fits of the VLT V-band light curve of OGLE-TR-56 (68 observations) from lPont et al.l 
ll2007ch for different treatment of LD. The parameter symbols are as given in Section[272] The transit midpoint To is expressed as HJD 
— 2 453 000. For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










TA + r b 


0.273 ±0.033 


0.275 ±0.035 


0.275 ±0.033 


0.275 ±0.031 


0.274 ±0.033 


k 


0.1022 ±0.0053 


0.1012 ±0.0056 


0.1015 ±0.0053 


0.1013 ±0.0047 


0.1017 ±0.0052 


i (deg.) 


78.8 ±2.5 


78.7 ±2.8 


78.7 ±2.7 


78.7±2.5 


78.8 ±2.6 


"A 


0.69 ±0.28 


0.49 ±0.33 


0.32 ±0.32 


0.79 ±0.35 


0.63 ±0.30 






0.25 perturbed 


0.55 perturbed 


0.23 perturbed 


0.10 perturbed 


T 


936.59816 ±0.00067 


936.59814 ±0.00066 


936.59813 ±0.00067 


936.59814 ±0.00064 


936.59814 ±0.00064 


TA 


0.248 ±0.030 


0.250 ±0.031 


0.250 ±0.030 


0.250 ±0.028 


0.249 ±0.030 


rb 


0.0253 ±0.0035 


0.0253 ±0.0037 


0.0254 ±0.0037 


0.0253 ±0.0034 


0.0253 ±0.0035 


cr (mmag) 


0.8418 


0.8429 


0.8428 


0.8429 


0.8425 


Y 2 

A red 


0.8062 


0.8083 


0.8080 


0.8083 


0.8075 


Fitting for two LD coefficients 


^A + r h 




0.258 ±0.015 


0.266 ±0.018 


0.259 ±0.015 


0.259 ±0.015 


k 




0.121 ±0.014 


0.113±0.010 


0.118±0.013 


0.118±0.012 


i (deg.) 




79.7 ±0.9 


79.1 ±1.1 


79.5 ±0.9 


79.5 ±0.9 


MA 




2.56 ±0.67 


3.96 ±1.80 


-0.26 ±0.37 


1.72 ±0.51 


^A 




-2.28 ±0.76 


-4.83 ±2.58 


-2.26 ±0.92 


-1.70 ±0.77 


To 




936.59827 ±0.00058 


936.59817 ±0.00062 


936.59826 ± 0.00056 


936.59827 ±0.00061 


TA 




0.230 ±0.014 


0.239 ±0.017 


0.232 ±0.015 


0.232 ±0.014 


rb 




0.0279 ±0.0028 


0.0270 ±0.0026 


0.0274 ±0.0026 


0.0273 ±0.0027 


cr (mmag) 




0.8229 


0.8241 


0.8219 


0.8224 


Y 2 

A rod 




0.7834 


0.7857 


0.7816 


0.7825 



Table A13. Parameters of the JKTEBOP best fits of the VLT i?-band light curve of OGLE-TR-56 (67 observations) from lPont et al.l 
(2007c) for different treatment of LD. The parameter symbols are as given in Section ^, 21 The transit midpoint To is expressed as HJD 
— 2 453 000. For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










^A + r h 


0.236 ±0.043 


0.265 ±0.039 


0.252 ±0.041 


0.347 ±0.018 


0. 333 ±0.042 


k 


0.0890 ±0.0104 


0.0912 ±0.0090 


0.0900 ±0.0099 


0.0918 ±0.0113 


0.0999 ±0.0081 


i (deg.) 


83.3 ±4.3 


81.4±3.3 


82.3 ±3.5 


73.4 ±1.5 


77.3 ±2.5 


MA 


1.14 ±0.31 


1.06 ±0.30 


1.06 ±0.29 


-0.84 ±3.43 


1.36±0.11 


^A 




0.28 perturbed 


0.60 perturbed 


0.25 perturbed 


0.10 perturbed 


To 


936.59811 ±0.00049 


936.59811 ±0.00049 


936.59814 ±0.00049 


936.59852 ±0.00046 


936.59843 ±0.00046 


TA 


0.217 ±0.037 


0.243 ±0.034 


0.231 ±0.035 


0.318 ±0.019 


0.303 ±0.036 


rb 


0.0193 ±0.0059 


0.0221 ±0.0051 


0.0208 ±0.0052 


0.0292 ±0.0023 


0.0302 ±0.0058 


cr (mmag) 


0.8531 


0.8514 


0.8521 


0.8452 


0.8550 


red 


0.8257 


0.8226 


0.8237 


0.9114 


0.8290 


Fitting for two LD coefficients 


rA + r h 




0.255 ±0.037 


0.225 ±0.037 


0.259 ±0.034 


0.257 ±0.034 


k 




0.0900 ±0.0099 


0.0910 ±0.0091 


0.0906 ±0.0092 


0.0904 ± 0.0096 


i (deg.) 




82.1 ±3.0 


83.7 ±2.4 


81.7 ±2.9 


81.9±2.9 


M A 




1.06 ±0.64 


1.54 ±1.42 


1.32 ±0.71 


1.13 ±0.38 


UA 




0.25 ±0.99 


-0.84 ±2.44 


0.37±1.02 


0.24 ±0.78 


T 




936.59811 ±0.00049 


936.59806 ±0.00050 


936.59807 ±0.00048 


936.59808 ±0.00048 


rA 




0.234 ±0.033 


0.206 ±0.033 


0.238 ±0.029 


0.236 ±0.030 


rb 




0.0210 ±0.0050 


0.0188 ±0.0049 


0.0215 ±0.0045 


0.0213 ±0.0045 


cr (mmag) 




0.8509 


0.8531 


0.8506 


0.8509 


X red 




0.8353 


0.8394 


0.8348 


0.8354 
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Table A14. Parameters of the JKTEBOP best fits of the OGLE-TR-111 light curve from IWinn et alJ j2007bh (I band, 7V obs = 386) for 
different treatment of LD. The parameter symbols are as given in Section[2]2] The transit midpoint Tq is expressed as HJD — 2 453 000. 
For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










TA + r b 


0.0954 ±0.0036 


0.0950 ±0.0040 


0.0954 ±0.0037 


0.0953 ±0.0036 


0.0953 ±0.0039 


k 


0.1324 ±0.0017 


0.1313 ±0.0021 


0.1318 ±0.0018 


0.1316 ±0.0018 


0.1314 ±0.0019 


i (deg.) 


88.03 ±0.47 


88.15 ±0.61 


88.08 ±0.51 


88.11±0.51 


88.12 ±0.56 


"A 


0.454 ±0.051 


0.320 ±0.064 


0.189 ±0.057 


0.604 ±0.061 


0.405 ± 0.060 






0.26 perturbed 


0.45 perturbed 


0.22 perturbed 


0.20 perturbed 


T 


799.75165 ±0.00019 


799.75165 ±0.00021 


799.75165 ±0.00020 


799.75165 ±0.00020 


799.75165 ±0.00020 


TA 


0.0843 ±0.0030 


0.0840 ±0.0034 


0.0843 ±0.0031 


0.0842 ±0.0031 


0.0842 ±0.0034 


rb 


0.01116 ±0.00053 


0.01103 ±0.00061 


0.01111 ±0.00053 


0.01108 ±0.00054 


0.01107 ±0.00058 


a (mmag) 


2.0830 


2.0828 


2.0828 


2.0827 


2.0828 


Y 2 

A red 


0.9786 


0.9778 


0.9781 


0.9779 


0.9779 


Fitting for two LD coefficients 


^A + r h 




0.0950 ±0.0040 


0.0954 ±0.0042 


0.0951 ±0.0040 


0.0951 ±0.0039 


k 




0.1313 ±0.0033 


0.1314 ±0.0040 


0.1312 ±0.0040 


0.1313 ±0.0047 


i (deg.) 




88.16 ±0.73 


88.11 ±0.63 


88.15 ±0.67 


88.15 ±0.64 


MA 




0.33 ±0.22 


-0.03 ±1.23 


0.66 ±0.51 


0.40 ±0.12 


^A 




0.24 ±0.44 


0.83 ±2.14 


0.29 ±0.70 


0.22 ±0.65 


To 




799.75165 ±0.00021 


799.75165 ±0.00020 


799.75165 ±0.00020 


799.75165 ±0.00019 


TA 




0.0839 ±0.0034 


0.0843 ±0.0036 


0.0841 ±0.0033 


0.0841 ±0.0033 


rb 




0.01102 ±0.00068 


0.01108 ±0.00066 


0.01103 ±0.00066 


0.01104 ±0.00064 


cr (mmag) 




2.0828 


2.0828 


2.0828 


2.0828 


Y 2 

A rod 




0.9804 


0.9806 


0.9805 


0.9805 



Table A15. Parameters of the JKTEBOP best fits of the OGLE-TR-132 light curve from lGillon et alJ ll2007ch (Special filter, N ohs = 274) 
for different treatment of LD. The parameter symbols are as given in Section ^. 21 The transit midpoint To is expressed as HJD — 2 453 000. 
For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










^A + r b 


0.228 ±0.017 


0.231 ±0.020 


0.231 ±0.018 


0.231 ±0.019 


0.231 ±0.020 


k 


0.0941 ±0.0021 


0.0937 ±0.0022 


0.0938 ±0.0020 


0.0937 ±0.0022 


0.0937 ±0.0023 


i (deg.) 


83.5 ±1.7 


83.2 ±2.0 


83.3 ±1.8 


83.3 ±2.0 


83.3 ±2.0 


MA 


0.534 ±0.071 


0.341 ±0.099 


0.166 ±0.088 


0.684 ±0.086 


0.464 ±0.088 


^A 




0.30 perturbed 


0.60 perturbed 


0.25 perturbed 


0.20 perturbed 


To 


142.59118 ±0.00033 


142.59121 ±0.00035 


142.59120 ±0.00033 


142.59121 ±0.00034 


142.59121 ±0.00034 


TA 


0.208 ±0.016 


0.212 ±0.017 


0.211 ±0.016 


0.211 ±0.017 


0.211 ±0.018 


rb 


0.0196 ±0.0018 


0.0198 ±0.0021 


0.0198 ±0.0019 


0.0198 ±0.0020 


0.0197 ±0.0021 


cr (mmag) 


1.2781 


1.2799 


1.2790 


1.2793 


1.2794 


Fitting for two LD coefficients 


r'A + r b 




0.226 ±0.016 


0.224 ±0.016 


0.225 ±0.017 


0.227±0.017 


k 




0.0944 ±0.0034 


0.0946 ±0.0034 


0.0944 ± 0.0034 


0.095 ±0.0037 


i (deg.) 




83.6 ±1.9 


83.7 ±2.1 


83.7±1.9 


83.5±1.8 


M A 




0.68 ±0.43 


1.13 ±1.28 


0.39 ±0.45 


0.61 ±0.22 


^A 




-0.22 ±0.65 


-0.96 ±2.09 


-0.25 ±0.76 


-0.19 ±0.59 


To 




142.59116 ±0.00034 


142.59115 ±0.00032 


142.59116 ±0.00033 


142.59118 ±0.00032 


TA 




0.206 ±0.014 


0.205 ±0.015 


0.206 ±0.015 


0.208 ±0.015 


rb 




0.0195 ±0.0019 


0.0194 ±0.0019 


0.0194 ±0.0020 


0.0196 ±0.0019 


cr (mmag) 




1.2801 


1.2803 


1.2801 


1.2800 
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Table A16. Parameters of the JKTEBOP best fits of the Spitzer light curve of GJ 436 for different LD laws. The parameter symbols are 
as given in Section[2]2] The transit midpoint To is expressed as HJD — 2454000. For each part of the table the upper parameters are 
fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


r'A + r b 


0.0794 ±0.0024 


0.0790 ±0.0024 


0.0794 ±0.0023 


0.0790 ±0.0023 


0.0791 ±0.0025 


k 


0.08313 ±0.00072 


0.08283 ±0.00075 


0.08259 ±0.00078 


0.08300 ±0.00077 


0.08291 ±0.00074 


i (deg.) 


86.41 ±0.15 


86.41 ±0.15 


86.41 ±0.15 


86.44 ±0.14 


86.44 ±0.15 


MA 


0.10 perturbed 


—0.05 perturbed 


—0.11 perturbed 


0.18 perturbed 


0.05 perturbed 






0.20 perturbed 


0.25 perturbed 


0.14 perturbed 


0.10 perturbed 


To 


280.78212 ±0.00011 


280.78211 ±0.00011 


280.78212 ±0.00011 


280.78212 ±0.00010 


280.78212 ±0.00010 




0.07332 ±0.0022 


0.0730 ±0.0022 


0.0734 ±0.0021 


0.0730 ±0.0021 


0.0730 ±0.0023 


rb 


0.006095 ±0.00019 


0.00604 ±0.00020 


0.00606 ±0.00018 


0.00606 ±0.00018 


0.00606 ±0.00019 


a (mmag) 


0.7043 


0.7051 


0.7041 


0.7049 


0.7048 


Y 2 

A red 


1.3718 


1.3748 


1.3712 


1.3743 


1.3738 
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Table A17. Parameters of the JKTEBOP best fits of the HD 149026 light curve from lSato et al.l l|2005h and lCharbonneau et al.l l|2006h for 
different treatment of LD. The parameter symbols are as given in Section ^, 21 The transit midpoint To is expressed as HJD — 2 453 000. 
For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 



Linear LD law Quadratic LD law Square-root LD law Logarithmic LD law 



APT (b+y/2) light curve from Sato et al. (2005) 



^A + r h 
k 


U.14J_ 019 

0464+ - 0033 
u.uwi_ 002g 


u.i4o_ 017 
0463+ 0033 

U.U4DJ_ 002g 


n 144+0-060 
U.144_ oi 7 

n n4«9+ - 0031 
0.0482_ 0035 


n 1 KO+0.065 
0-15»_ .015 

n n4sn+ 0033 
0.0480_ 0030 


i (deg.) 


82.0 to 90.0 


83.0 to 90.0 


83.0 to 90.0 


83.0 to 90.0 


«A 


yjer+O- 35 

u./oo_ 26 


U.OlU_ 302 


n +0-358 
u.ooi_ 271 


0.780l°; 29 « 


VA 




0.25 perturbed 


0.55 perturbed 


0.20 perturbed 


P (days) 
T 


9 S7W+0 00037 

527.87291°;^ 


2.87560lHog 3 36 
527.87291°™" 


9 07^8+0-00034 
_.OfOOO_ 00036 

527 8729+ 0014 
OZ '-°' ZM -0.0015 


9 S79SK+0-00031 

_.8/_8b_ a00034 
527.87291°;°™ 


TA 

rb 


0.136lHf 8 
00633+ 0039 

U.UUDOO_q Q009 


n ,00+0.062 
0064+ - 0009 

U.UUDt_ Q 0040 


iro+0.057 
U.lOO_ oi6 

0072+ ' 0037 
u.uu(z_ Q 00Q8 


o.i5ll°;°? 4 

0071+ 0041 
u.uun_ Q007 


a (mmag) 


1.8014 


1.8022 


1.8018 


1.8019 


FLWO g light curve from Charbonncau et al. (2006) 


^A + r h 
k 


162+ 045 

U.1DZ_ 014 

0491+ 0039 
U,U43 -0 .0027 


n ic;s+o.oio 
u.ioo_ 004 

n533+°- 0018 
u.uooo_ 0014 


0.158tooQ 4 
0.053812-°™ 


n 1 47+0.008 

u - ±z± < -0.004 

n n4s^+ - 0025 
0.0485_ 0024 


i (deg.) 


82.0 to 90.0 


86.0 to 90.0 


86.0 to 90.0 


87.0 to 90.0 


u A 


0.73 perturbed 


0.52 perturbed 


0.42 perturbed 


0.80 perturbed 


VA 




0.24 perturbed 


0.41 perturbed 


0.20 perturbed 


To 


527.8750+^i 3 


527.8749+™ 


527.8749lg°°i 3 


527.875ll° : °ol2 


rA 
rb 


0.154+ToiS 

n075 +0 - 0009 
U - UU '°-0.0030 


o 1 5n+ 0034 

u.xou_ oogg 
0.00801°;°°°? 


0-150l°°» 

n nnsn+ - 0009 
0.0080_ 0003 


0-i40l°;°g 4 

0-00681™ 


a (mmag) 


2.4386 


2.2049 


2.2054 


2.4354 


V 2 

A red 


1.0053 


1.0019 


1.0027 


1.0018 


FLWO r li 


ght curve from Charbonncau ct al. (2006) 






r A + r b 
k 


162+ 045 

U.1DZ_ Q14 

0491+ - 0039 
-0.0027 


147+0.009 
u - 1 ^' -0.004 

0484+ - 0023 
u - u ^°^-0.0023 


0-148_ 004 
0.0486_ Q 0024 


147 +0.008 
-0.004 

n n4s=;+ - 0025 
0.0485_ 0024 


i (deg.) 


82.0 to 90.0 


87.0 to 90.0 


85.0 to 90.0 


87.0 to 90.0 


u A 


0.61 perturbed 


0.35 perturbed 


0.19 perturbed 


0.70 perturbed 


VA 




0.32 perturbed 


0.55 perturbed 


0.27 perturbed 


To 


527.8750t°- 99 » 


527.8750+H 9 ! 2 


527.875ll°;~i 2 2 


527.875ll°;»i2 


rA 
rb 


0.154l° ;°g 
0076+ - 0028 

U.UU(D_ 00Q9 


0-140l°;°°« 
n nnfis+ - 0006 


0.1411°;^ 


0.1401°;°°* 
0.0068_ 0004 


a (mmag) 


2.4386 


2.4362 


2.4353 


2.4354 


^ red 


1.0053 


1.0027 


1.0017 


1.0018 
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Table A18. Parameters of the JKTEBOP best fits of the FLWO r-band light curve (jV oba = 642) of HD 189733 from lBakos et al.l (|2006bl l 
for different treatments of LD. The parameter symbols are as given in Section[2]2] The transit midpoint To is expressed as HJD — 
2 453 000. For each part of the tabic the upper parameters are fitted quantities and the lower parameters are derived quantities. Note 
that the a is affected by scatter due to high airmass towards the end of the observations. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










ta + r h 


0.1190 ±0.0068 


0.1194 ±0.0066 


0.1193 ±0.0069 


0.1194 ±0.0069 


0.1192 ±0.0073 


k 


0.1509 ±0.0034 


0.1500 ±0.0036 


0.1504 ±0.0037 


0.1502 ±0.0033 


0.1497 ±0.0039 


i (deg.) 


86.56 ±0.61 


86.56 ±0.64 


86.55 ±0.65 


86.55 ±0.64 


86.59 ±0.69 




0.541 ±0.123 


0.413 ±0.148 


0.309 ±0.133 


0.651 ±0.139 


0.481 ±0.148 


VA 




0.21 perturbed 


0.38 perturbed 


0.18 perturbed 


0.20 perturbed 


To 


642.70584 ±0.00025 


642.70584 ±0.00026 


642.70584 ±0.00025 


642.70584 ±0.00024 


642.70584 ±0.00025 


TA 


0.1034 ±0.0056 


0.1038 ±0.0055 


0.1037 ±0.0057 


0.1038 ±0.0057 


0.1037 ±0.0060 




0.0156 ±0.0011 


0.0156 ±0.0012 


0.0156 ±0.0012 


0.0156 ±0.0012 


0.0155 ±0.0013 


cr (mmag) 


4.9239 


4.9249 


4.9243 


4.9245 


4.9249 


Fitting for two LD coefficients 


+ r h 




0.1184 ±0.0052 


0.1181 ±0.0059 


0.1182 ±0.0064 


0.1183 ±0.0059 


k 




0.1542 ±0.0079 


0.1536 ±0.0072 


0.1539 ±0.0093 


0.1537 ±0.0080 


i (deg.) 




86.50 ±0.50 


86.52 ±0.52 


86.51 ±0.58 


86.50 ±0.55 


MA 




0.98 ±0.75 


1.71 ±2.23 


0.11 ±1.01 


0.73 ±0.35 


^A 




-0.70 ±0.49 


-1.90±3.63 


-0.73 ±1.50 


-0.53 ±1.07 


To 




642.70584 ±0.00025 


642.70584 ±0.00025 


642.70584 ±0.00025 


642.70584 ±0.00025 


TA 




0.1026 ±0.0042 


0.1024 ±0.0050 


0.1024 ±0.0054 


0.1026 ±0.0049 


rh 




0.0158 ±0.0011 


0.0157 ±0.0011 


0.0158 ±0.0012 


0.0158 ±0.0012 


cr (mmag) 




4.9227 


4.9231 


4.9229 


4.9230 



Table A19. Parameters of the JKTEBOP best fits of the FLWO light curve of HD 189733 (z band, Af obs = 1662) from lWinn et al.l d2007dh 
for different treatments of LD. The parameter symbols are as given in Section ^, 21 For each part of the table the upper parameters are 
fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










^a + n> 


0.1295 ±0.0019 


0.1298 ±0.0021 


0.1297 ±0.0019 


0.1298 ±0.0020 


0.1297 ±0.0020 


k 


0.1585 ±0.0009 


0.1572 ±0.0009 


0.1579 ±0.0009 


0.1576 ±0.0009 


0.1580 ±0.0009 


i (deg.) 


85.73 ±0.16 


85.73 ±0.16 


85.73 ±0.16 


85.74 ±0.17 


85.73 ±0.17 


M A 


0.440 ±0.058 


0.222 ±0.062 


0.124 ±0.062 


0.571 ±0.066 


0.397 ±0.058 


«A 




0.31 perturbed 


0.50 perturbed 


0.24 perturbed 


0.10 perturbed 


^A 


0.1118 ±0.0016 


0.1122 ±0.0016 


0.1120 ±0.0016 


0.1121 ±0.0017 


0.1120 ±0.0016 


rh 


0.01772 ±0.00033 


0.01764 ±0.00033 


0.01768 ±0.00033 


0.01767 ±0.00035 


0.01769 ±0.00034 


a (mmag) 


2.4579 


2.4578 


2.4578 


2.4578 


2.4578 


^ red 


1.0086 


1.0085 


1.0086 


1.0085 


1.0086 


Fitting for two LD coefficients 


^a + n> 




0.1299 ±0.0023 


0.1298 ±0.0026 


0.1299 ±0.0027 


0.1297 ±0.0026 


k 




0.1573 ±0.0043 


0.1575 ±0.0044 


0.1575 ±0.0050 


0.1578 ±0.0048 


i (deg.) 




85.73 ±0.18 


85.74 ±0.17 


85.73 ±0.19 


85.73 ±0.18 


M A 




0.23 ±0.65 


-0.02 ±1.79 


0.59 ±0.72 


0.37±0.29 


VA 




0.30 ±0.97 


0.73 ±2.91 


0.28 ±1.20 


0.15 ±0.84 


VA 




0.1122 ±0.0022 


0.1121 ±0.0024 


0.1122 ±0.0026 


0.1121 ±0.0025 


rh 




0.01765 ±0.00045 


0.01777 ±0.00038 


0.01767 ±0.00041 


0.01768 ±0.00039 


a (mmag) 




2.4578 


2.4578 


2.4578 


2.4578 


^ red 




1.0091 


1.0092 


1.0091 


1.0091 
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Table A20. Parameters of the JKTEBOP best fits of the WISE Observatory HD 189733 light curve (I band, 7V obs = 345) from I Winn et al.l 
[|2007dh for different LD laws. The parameter symbols are as given in Section[2]2] The transit midpoint To is expressed as HJD — 
2 450 000. For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law Lc 


igarithmic LD law 


Cubic LD law 


Fitting for one LD coefficient 


r\ + r b 


0.1320 ±0.0067 


0.1327 ±0.0072 


0.1327 ±0.0073 


0.1327 ±0.0071 


0.1325 ±0.0074 


k 


0.1506 ±0.0025 


0.1496 ±0.0026 


0.1500 ±0.0024 


0.1498 ±0.0025 


01501 ±0.0026 


i (deg.) 


85.41 ±0.56 


85.37 ±0.61 


85.38 ±0.58 


85.38 ±0.60 


85.38 ±0.60 


"A 


0.42 ±0.22 


0.22 ±0.27 


0.11±0.25 


0.52 ±0.27 


0.36 ±0.25 






0.25 perturbed 


0.46 perturbed 


0.22 perturbed 


0.10 perturbed 


T 


984.36572 ±0.00029 


984.36582 ±0.00028 


984.36582 ±0.00027 984.36582 ± 0.00028 984.36582 ± 0.00029 


TA 


0.1147 ±0.0056 


0.1154 ±0.0060 


0.1153 ±0.0062 


0.1154 ±0.0060 


0.1152 ±0.0064 




0.0173 ±0.0010 


0.0173 ±0.0011 


0.0173 ±0.0011 


0.0173 ±0.0011 


0.0173 ±0.0011 


a (mmag) 3.2212 


3.2204 


3.2207 


3.2205 


3.2207 


Y 2 

A red 


1.0042 


1.0038 


1.0039 


1.0038 


1.0040 


Table A21. 


Parameters of the JKTEBOP best fits of the HD 209458 light curve from |b 


rown et alJ J200lf) (N„u. = 5561 for different 


treatments of LD. The parameter symbols are as given in Section^, 21 The transit midpoint Tq is expressed as HJD — 2 451 000. For each 


part of the table the upper parameters 


are fitted quantities and the lower parameters are 


derived quantities. 






Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










rx + r h 


0.12889 ±0.00041 


0.12768 ±0.00048 


0.12840 ±0.00044 


0.12811 ±0.00048 


0.12846 ±0.00043 


k 


0.12260 ±0.00011 


0.12095 ±0.00020 


0.12172 ±0.00013 


0.12137 ±0.00018 


0.12201 ±0.00021 


i (deg.) 


86.472 ±0.038 


86.668 ±0.052 


86.567 ±0.044 


86.612 ±0.049 


86.543 ±0.046 


MA 


0.494 ±0.004 


0.297 ±0.020 


0.152 ±0.021 


0.674 ±0.023 


0.463 ±0.011 


VA 




0.34 perturbed 


0.57 perturbed 


0.28 perturbed 


0.10 perturbed 


To 


659.936716 ± 0.000021 


659.936712 ± 0.000022 


659.936714 ±0.000021 


659.936713 ±0.000020 


659.963715 ±0.000022 


VA 


0.11482 ±0.00035 


0.11391 ±0.00042 


0.11447 ±0.00038 


0.11425 ±0.00041 


0.11449 ±0.00038 


rb 


0.014076 ± 0.000055 


0.013778 ± 0.000068 


0.013934 ±0.000059 


0.013866 ±0.000065 


0.013969 ±0.000062 


a (mmag) 


0.1781 


0.1732 


0.1748 


0.1738 


0.1754 


Y 2 
red 


1.1457 


1.0536 


1.0812 


1.0633 


1.0947 


Fitting for two LD coefficients 


rA + r b 




0.12769 ±0.00047 


0.12794 ±0.00050 


0.12779 ±0.00053 


0.12774 ±0.00050 


k 




0.12096 ±0.00025 


0.12052 ±0.00034 


0.12073 ±0.00030 


0.12074 ±0.00032 


i (deg.) 




86.667 ±0.53 


86.688 ±0.057 


86.681 ±0.063 


86.689 ±0.060 


MA 




0.297 ±0.028 


-0.309 ±0.122 


0.772 ±0.042 


0.405 ±0.013 


VA 




0.338 ±0.050 


1.350 ±0.207 


0.427 ±0.063 


0.317±0.050 


To 




659.936712 ± 0.000022 


659.93671 ±0.000022 


659.936711 ±0.000023 


659.936711 ±0.000022 


rA 




0.11391 ±0.00041 


0.11418 ±0.00043 


0.11402 ±0.00045 


0.11398 ±0.00043 


rb 




0.013779 ± 0.000069 


0.013761 ±0.000075 


0.013766 ±0.000075 


0.013762 ±0.000078 


a (mmag) 




0.1732 


0.1734 


0.1733 


0.1734 


^ red 




1.0536 


1.0540 


1.0538 


1.0541 
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Table A22. Parameters of the JKTEBOP best fits of the MOST light curve of HD 209458 from iRowe et alj |2006) (unfiltered GOAT?, 
N a b s = 8772 as analysed here) for different treatments of LD. The parameter symbols are as given in Section ^, 21 For each part of the 
table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


Fitting for 


one LD coefficient 










r\ + r b 


0.1269 ±0.0021 


0.1259 ±0.0022 


0.1266 ±0.0022 


0.1264 ±0.0022 


0.1266 ±0.0022 


k 


0.12447 ±0.00071 


0.1231 ±0.00079 


0.1235 ±0.00080 


0.1235 ±0.00078 


0.1238 ±0.00076 


i (deg.) 


86.69 ±0.22 


86.85 ±0.24 


86.77 ±0.24 


86.79 ±0.24 


86.75 ±0.23 


u A 


0.552 ±0.023 


0.415 ±0.031 


0.196 ±0.032 


0.691 ±0.034 


0.524 ±0.026 


VA 




0.25 perturbed 


0.60 perturbed 


0.21 perturbed 


0.10 perturbed 


VA 


0.1128 ±0.0018 


0.1121 ±0.0019 


0.1127 ±0.0019 


0.1125 ±0.0019 


0.1126 ±0.0019 


rb 


0.01404 ±0.00030 


0.01380 ±0.00031 


0.01392 ±0.00032 


0.01389 ±0.00031 


0.01395 ±0.00030 


a (mmag) 


3.1582 


3.1564 


3.1571 


3.1570 


3.1574 


% red 


0.9258 


0.9247 


0.9251 


0.9251 


0.9253 


Fitting for two LD coefficients 


+ n> 




0.1292 ±0.0058 


0.1312 ±0.0048 


0.1321 ±0.0055 


0.1321 ±0.0055 


k 




0.1122 ±0.0014 


0.1162 ±0.0020 


0.1142 ±0.0012 


0.1151 ±0.0023 


i (deg.) 




87.46 ±0.56 


86.95 ±0.36 


87.02 ±0.43 


86.95 ±0.39 


u A 




-0.05 ±0.19 


-2.15±0.73 


2.08 ±0.26 


0.38 ±0.07 


v A 




1.724 ±0.36 


4.82 ±1.31 


1.94 ±0.38 


1.39 ±0.39 


rA 




0.1162 ±0.0053 


0.1175 ±0.0045 


0.1185 ±0.0050 


0.1185 ±0.0052 


rb 




0.01304 ±0.00055 


0.01366 ±0.00045 


0.01354 ± 0.00050 


0.01364 ±0.00048 


a (mmag) 




3.1524 


3.1542 


3.1534 


3.1538 


^ red 




0.9228 


0.9238 


0.9234 


0.9236 



Table A23. Parameters of the JKTEBOP best fits of the HST light curves of HD 209458 from iKnutson et alj <2007d) (JV obs = 504 for 
the five bluer passbands and -/V ODS = 548 for the five redder passbands) for the linear LD law. The parameter symbols are as given in 
Section ^. 21 For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 





320 nm 


375 nm 


430 nm 


484 nm 


539 nm 


»"A + r b 


0.1270 ±0.0024 


0.1247 ±0.0024 


0.1271 ±0.0007 


0.1272 ±0.0006 


0.1286 ±0.0006 


k 


0.12230 ±0.00128 


0.12160 ±0.0013 


0.12248 ±0.0003 


0.12238 ±0.0002 


0.12272 ±0.0002 


i (deg.) 


86.71 ±0.26 


86.93 ±0.15 


86.65 ±0.07 


86.63 ±0.06 


86.50 ±0.06 


u A 


0.828 ±0.023 


0.754 ±0.013 


0.703 ±0.007 


0.618 ±0.006 


0.561 ±0.0007 


rA 


0.1131 ±0.0021 


0.1112±0.0011 


0.1133 ±0.0006 


0.1134 ±0.0005 


0.1146 ±0.0005 


rb 


0.01384 ±0.00039 


0.01352 ±0.00020 


0.01387 ±0.00010 


0.01387 ±0.00008 


0.01406 ±0.00009 


a (mmag) 


0.6568 


0.3686 


0.2238 


0.1989 


0.2067 


^ red 


1.0791 


1.0700 


1.0129 


1.0843 


1.0264 




580 nm 


677 nm 


775 nm 


873 nm 


970 nm 


r\ + r b 


0.1286 ±0.0005 


0.1296 ±0.0004 


0.1281 ±0.0006 


0.1294 ±0.0008 


0.1290 ±0.0010 


k 


0.12328 ±0.00017 


0.12266 ±0.00013 


0.12222 ±0.00015 


0.12270 ±0.00042 


0.12336 ±0.00024 


i (deg.) 


86.48 ±0.05 


86.38 ±0.04 


86.49 ±0.06 


86.37 ±0.07 


86.39 ±0.10 


UA 


0.534 ±0.006 


0.437 ±0.006 


0.377 ±0.008 


0.324 ±0.011 


0.275 ±0.016 


rA 


0.1145 ±0.0004 


0.1154 ±0.0004 


0.1141 ±0.0005 


0.1152 ±0.0007 


0.1148 ±0.0009 


rb 


0.01411 ±0.00007 


0.01416 ±0.00006 


0.01395 ± 0.00008 


0.01413 ±0.00010 


0.01416 ±0.00013 


a (mmag) 


0.2210 


0.1861 


0.2413 


0.3161 


0.4912 


^ red 


1.3139 


1.0623 


1.0866 


1.0341 


1.1273 
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Table A24. Parameters of the JKTEBOP best fits of the HST light curves of HD 209458 from lKnutson et al.l [|2007ch (N obs = 504 for the 
five bluer passbands and N obs = 548 for the five redder passbands) for the quadratic LD law. The parameter symbols are as given in 
Section ^. 21 For each part of the table the upper parameters are fitted quantities and the lower parameters are derived quantities. 







320 nm 


375 nm 


430 nm 


484 nm 


539 nm 


r-A + r b 
k 

i (deg.) 

UA 
VA 




0.1271 ±0.0022 
0.1249 ±0.0018 
86.57 ±0.25 
1.030 ±0.102 
-0.384 ±0.182 


0.1248 ±0.0012 
0.1221 ±0.0009 
86.90 ±0.15 
0.791 ±0.052 
-0.073 ±0.102 


0.1271 ±0.0007 
0.12249 ±0.0006 
86.65 ±0.09 
0.703 ±0.036 
-0.001 ±0.068 


0.1272 ±0.0006 
0.1223 ±0.0005 
86.64 ±0.07 
0.612 ±0.034 
0.009 ±0.062 


0.1277 ±0.0008 
0.1213 ±0.0005 
86.66 ±0.09 
0.426 ±0.039 
0.248 ±0.092 


TA 

rb 

a (mma; 

red 





0.1129 ±0.0018 
0.01411 ±0.0004 
0.6528 
1.0695 


0.1112±0.0010 
0.01357 ±0.00021 
0.3683 
1.0704 


0.1133 ±0.0006 
0.01387 ±0.00013 
0.2238 
1.1315 


0.1133 ±0.0005 
0.01386 ±0.00010 
0.1989 
1.0862 


0.1139 ±0.0007 
0.01382 ±0.00013 
0.2038 
1.0027 






580 nm 


677 nm 


775 nm 


873 nm 


970 nm 


ta + r h 
k 

i (deg.) 

UA 
VA 




0.1281 ±0.0006 
0.1227 ±0.0004 
86.55 ±0.07 
0.462 ±0.036 
0.126 ±0.063 


0.1288 ±0.0005 
0.1218 ±0.0003 
86.49 ± 0.06 
0.309 ±0.037 
0.214 ±0.061 


0.1271 ±0.0006 
0.1211 ±0.0003 
86.64 ±0.07 
0.197 ±0.047 
0.299 ±0.078 


0.1283 ±0.0008 
0.1214 ±0.0004 
86.54 ±0.09 
0.079 ±0.069 
0.400 ±0.114 


0.1272 ±0.0013 
0.1215 ±0.0006 

86.65 ±0.16 
-0.078 ±0.098 

0.581 ±0.164 


rA 
rb 

a (mma; 

^ red 





0.1141 ±0.0005 
0.01400 ±0.00009 
0.2206 
1.3088 


0.1148 ±0.0004 
0.0140 ±0.00008 

0.1846 

1.0429 


0.1134 ±0.0005 
0.01373 ±0.00009 
0.2382 
1.0593 


0.1144 ±0.0007 
0.01388 ±0.00012 
0.3122 
1.0105 


0.1135 ±0.0011 
0.01379 ±0.00019 
0.4861 
1.106 
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